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1.0  EXECUTIVE  SUMMARY 


This  report  summarizes  the  experimental  and  theoretical  research  effort  undertaken  to  study  and 
identify  some  of  the  chemical  pathways  involved  in  the  hypergolic  ignition  of  mixtures 
consisting  of  room  temperature  ionic  liquids  (ILs)  and  suitable  oxidizers  such  as  nitric  acid. 
Hypergolicity  is  defined  as  the  auto-ignition  of  a  binary  chemical  system  under  ambient 
conditions  of  pressure  and  temperature.  The  research  was  motivated  by  the  recognition  that 
higher  performing  hypergolic  propellants  need  to  be  developed  to  meet  future  DoD  space 
propulsion  capabilities  as  outlined  in  the  IHPRPT  (Integrated  High  Payoff  Rocket  Propulsion 
Technology)  Program.  The  popularity  of  ILs  as  a  class  of  energetic  substances  for  DoD 
applications  arises  from  the  fact  that  a  combination  of  different  anions  and  cations  in  the  IL  will 
lead  to  a  large  (-10s)  number  of  compounds  with  tailorable  physical  and  thermochemical 
properties.  While  many  physical  and  thermodynamic  properties  of  propellants  are  now  being 
predicted  theoretically  with  reasonable  reliability  using  quantum  ab  initio  computations  and 
molecular  dynamics  simulations,  the  same  is  not  true  for  a  priori  predictions  of  hypergolic 
reactivity.  Our  aim  in  this  work  was  to  provide  accurate  mechanistic  data  on  recently  discovered 
hypergolic  ILs  that  can  be  used  to  build  a  chemical  kinetics  ignition  model  to  simulate  the 
ignition  delay  time,  a  quantity  of  significant  importance  in  the  design  and  implementation  of 
hypergolic  rocket  engines.  The  mechanistic  information  can  also  be  used  for  developing  reactive 
molecular  dynamics  models  and  correlation  algorithms  to  create  “first-principle”  design  tools  to 
predict  new  hypergolic  combinations  prior  to  any  effort  to  synthesize  the  fuel  molecules. 
Availability  of  such  computational  tools  would  provide  a  mitigation  strategy  for  the  risks  (cost 
and  time)  currently  involved  in  the  trial-and-error  approach  to  first  synthesize  the  IL  and  then 
test  it  to  verify  its  behavior.  An  added  feature  of  hypergolic  IL  fuels  is  their  inherent  low  vapor 
pressures.  Implementation  of  such  fuels  would  result  in  lowering  of  vapor  toxicity  exposures 
and  reduction  in  emissions  to  the  environment.  Consequently,  fielding  of  such  “green”  IL  fuels 
in  future  propulsion  systems  will  provide  significant  logistical,  operational  and  supportability 
cost  savings  to  the  DoD. 

The  details  of  our  findings  can  be  found  in  Appendix  A1'4  and  in  various  publications5'8  as  well 
as  in  symposia  proceedings.9'11  Below  is  a  brief  synopsis  of  ongoing  unpublished  work. 


2.0  THEORETICAL  EFFORT 
2.1  Direct  Dynamics  Trajectory  Studies 

Based  on  experimental  investigations  of  the  end  product  analysis,  a  possible  reaction  mechanism, 
involving  the  1,5-dinitrobiuret  (DNB)  intermediate,  in  the  hypergolic  ignition  of  dicyanamide- 
based  ILs  and  nitric  acid  system  was  published.1  To  further  understand  the  chemistry  of  this 
nitrogen-rich  explosive  compound,  ab  initio  electronic  structure  calculations  and  direct  dynamics 
simulations  were  carried  out  to  map  the  reaction  coordinates  for  decomposition  of  DNB, 
including  energies  and  structures  of  all  stationary  points  (i.e.,  intermediate  complexes,  transition 
states  (TSs),  and  products).  For  complex  molecules,  especially  for  those  with  high  internal 
energy  content,  the  use  of  chemical  intuition  to  predict  reaction  paths  is  unreliable.  There  may 
be  concerted  reactions  that  are  difficult  to  predict,  and  once  energy  stored  in  the  molecule  begins 
to  release,  the  system  does  not  necessarily  follow  the  minimum  energy  reaction  path,  i.e.,  the 
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subsequent  behavior  is  controlled  by  dynamics.  Even  for  steps  in  the  reaction  that  are  controlled 
by  statistical  factors,  standard  transition  state  locating  methods  may  be  difficult  to  apply,  because 
the  molecular  transformations  are  too  complex.  A  useful  approach  to  treating  such  a  system  is 
quasi-classical,  direct  dynamics  trajectory  simulations,  where  the  motion  of  the  molecule  is 
followed,  allowing  the  molecule  to  "show  us"  what  the  preferred  reaction  paths  are.  The  direct 
dynamics  method  dispenses  with  the  potential  energy  surface.  Instead,  it  calculates  the  energies, 
force  constants,  and  Hessian  "on  the  fly"  using  ab  initio  methods.  This  method  becomes 
computationally  attractive  when  the  dimensionality  of  the  system  increases,  particularly  for  the 
DNB  molecule  which  contains  13  heavy  atoms.  Dynamics  simulations  partition  the  energy 
generated  by  the  exothermic  reaction  to  vibrational,  rotational  and  translational  modes, 
increasing  the  chance  of  locating  new  reaction  paths.  In  addition,  by  following  the  variation  of 
the  potential  energy  during  the  trajectory  rather  than  relying  on  intuition,  we  can  identify  better 
geometries  for  TS  searching. 

In  this  study,  direct  dynamics  trajectory  simulations  were  conducted  for  two  purposes:  1) 
examining  decomposition  behavior  of  DNB  at  various  temperatures  and  discovering  new 
decomposition  paths;  and  2)  locating  key  transition  states  for  the  dissociation  paths  identified  in 
dynamics  simulations,  when  the  standard  TS-searching  methods  (i.e.,  TS,  QST2,  and  QST3 
methods  in  Gaussian  03 12)  fail.  Trajectory  simulations  were  started  at  the  equilibrium  geometry 
of  DNB.  A  general  chemical  dynamics  program  VENUS13  was  used  to  set  up  the  trajectory 
initial  conditions,  and  the  Bakken  Hessian-based  method  implemented  in  Gaussian  03  was  used 
to  propagate  each  trajectory  with  Hessians  recalculated  every  five  steps.  The  integrators  were 
performed  with  a  step  size  of  ~  0.4  fsec,  which  conserved  total  energy  to  better  than  10'4  Hartree. 
The  scf  =  xqc  option  was  used  during  trajectory  integration  so  that  a  quadratically  convergent 
Hartree-Fock  (qc-scf)  method  was  used  in  case  the  usual,  but  much  faster  scf  method  did  not 
converge  within  the  allotted  number  of  cycles.  B3LYP/6-3 1G*  level  of  theory  was  chosen  for 
the  main  set  of  trajectories.  For  most  of  the  trajectories  calculated  in  the  vibrational  temperature 
range  4000-6000  K,  the  time  at  which  the  DNB  began  to  dissociate  was  between  0.5  and  3  psec, 
which  is  much  longer  than  the  periods  of  most  DNB  vibrational  modes.  At  a  vibrational 
temperature  of  4000  K,  10%  of  trajectories  maintain  the  original  molecular  structure  before 
trajectory  termination.  The  dissociative  trajectories  can  be  grouped  into  six  classes,  as  illustrated 
in  Table  1. 
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Table  1.  Direct  Dynamics  Simulation  Results  for  Thermal  Decomposition  of  DNB 


Dissociation  Channels 


AH a  Barrier  of  Branching  % b 

/eV  possible  TS  /  eV  4000 K  I  5000K  I  6000 K 


Path  1: 


%N/°  °\N^°  V/° 

i  T  i  — ►  i  i 
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I  H  Ox  yO 

l  k  n  1.81  TS1, 0.63 


o  o 


o  o 


36  41  37 


▼  H 
°%l^°  N 

l  +  c"  x2 
/N  II 

H  O 


note:  some  trajectories  dissociate  simultaneously 
Path  2: 
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T  i  i 


V^° 


/u\  1.15 

H  1ST 


TS2,  1.22 


24  12 
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/°\  +  N  -0.59 
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Path  3: 
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N\  u/NVNVN\/° 
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ii  it  ll  ii  / 

U  o  o  O  o7 

°%-°  °-N^° 

I  +  I 

N ..  ,N, 

k'  V  X.  nh 


0.66  TS3, 0.36 
(TS3  corre¬ 
sponds  to 
twisting  CN 
bond  before  H- 
transfer) 
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Table  1.  continued. 


Dissociation  Channels 


Path  4: 


AH a  Barrier  of  Branching  % 

/eV  possible  TS  /  eV  4000K  I  5000K  I  . 


°s//°  %N/° 

I  'l'  l  — ►  l  i 


<W° 


,/N\0/N\ /N\ 


O  o 


h  V  V  'h  h  'C  c  H 


o  .  O 


X,/0 

I  I 

h'  VN-SN. 


H  N  0.97 


‘S>  0*  O 

V 

I 
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Path  5: 
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H^vV\  /vV\ 
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o  O 


o  .  O 


6  17 


<W0  +  # 


+  /\  +  c 

NH - NH  H 


Path  6: 


%  /°  °N  %  /O 

i  i  i  —  i  i  ♦  v°- 

H/Vc/Nx\  h/nxnxn^h 


%/° 

V  H  M 

I  |  II  +  NO 

./WS, 


Non-reactive  within  4  psec  simulation  period  0  10  - 

Sum  of  other  minor  dissociation  paths  -  6  15 

aValues  of  AH  were  calculated  relative  to  DNB  at  the  B3LYP/6-31G*  level  of  theory. 

branching  ratios  were  calculated  based  on  50,  60  and  30  trajectories  simulated  at  4000,  5000,  and  6000  K, 

respectively. 
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Ab  initio  electronic  structure  calculations  were  performed  at  the  B3LYP/6-31G*  level  of  theory, 
using  the  Gaussian  03  suite  of  programs.  Geometries  were  optimized  by  calculating  the  force 
constants  at  every  step.  Vibrational  frequencies  and  zero-point  energies  (ZPEs)  were  scaled  by  a 
factor  of  0.9613  and  0.9804,  respectively,  and  when  necessary  intrinsic  reaction  coordinate  (IRC) 
calculations  were  used  to  determine  which  minima  are  connected  by  particular  TSs.  All 
transition  states  found  were  verified  to  be  first-order  saddle  points  by  frequency  calculation. 
Optimized  geometries  of  DNB  and  various  dissociation  products  identified  in  the  direct 
dynamics  simulations  are  depicted  in  Fig.  1,  and  those  of  transition  states  for  the  three  major 
dissociation  paths  are  depicted  in  Fig.  2. 


(e) 

y 


Fig.  1:  B3LYP/6-31G*  optimized  geometries  of  the  reactant  and  products  identified  in 
trajectory  simulations:  a)  DNB,  b)  02N-NHC0-NHC0-NH,  c)  02N-NHC0-NHC0,  d)  02N- 
NHCO-NCO,  e)  02N-NHC0-NHC0-N0-H,  f)  cyclic  [HN-CO-NH],  g)  HNN(0H)0,  h) 
H2NN02,  i)  HNN02,  j)  HNCO,  k)  N02,  and  1)  N20. 
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Fig.  2:  B3LYP/6-31G*  optimized  geometries  of  transition  states.  Vibrational  modes 
corresponding  to  the  imaginary  frequencies  were  indicated  by  displacement  vectors. 


By  far  the  most  common  outcome  of  DNB  decomposition  is  Path  1,  where  elimination  of  NO2 
occurs  first,  followed  by  sequential  dissociation  of  the  remaining  unit  to  isocyanic  acid  HNCO 
and  HNNO2.  A  typical  Path  1  trajectory  (calculated  at  4000  K  vibrational  temperature)  is 
illustrated  in  Fig.  3,  showing  the  change  of  potential  energy  (PE)  during  the  trajectory.  The 
oscillations  in  the  PE  reflect  the  vibration  of  the  DNB  molecule,  including  ZPE.  The  trajectory 
shows  three  dissociation  steps,  eliminating  NO2  followed  by  elimination  of  two  HNCO 
molecules.  Elimination  of  NO2  occurs  at  ~  150  fsec,  and  the  time  for  which  the  decomposition 
lasts  is  ~  650  fsec. 
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Fig.  3:  A  representative  trajectory  of  dissociation  Path  1  at  4000  K,  simulated  at  B3LYP/6- 
31G*.  Plot  shows  the  variation  of  potential  energy  during  the  trajectory,  and  three 
dissociation  steps  are  indicated. 


In  Path  2,  HNN(OH)0  is  produced  via  an  intra-molecular  hydrogen  atom  transfer  followed  by 
dissociation  of  an  amide  bond.  HNN(OH)0  further  decomposes  to  N20  and  water.  A  trajectory 
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representative  of  dissociation  Path  2  is  illustrated  in  Fig.  4,  indicating  two  energy  barriers 
associated  with  two  activated  dissociation  steps. 


Trajectory  Time  /  fsec 


Fig.  4:  A  representative  trajectory  of  dissociation  Path  2  at  4000  K,  calculated  at  B3LYP/6- 
31G*.  Plot  shows  two  dissociation  steps  and  the  variation  of  potential  energy  during  the 
trajectory. 


The  other  common  dissociation  channel  is  Path  3  where  DNB  first  twists  along  a  CN  bond 
through  TS3  (see  Fig.  2),  then  undergoes  an  intra-molecular  hydrogen  atom  transfer  (i.e.,  transfer 
of  an  H-atom  from  the  central  nitrogen  atom  to  the  neighboring  N)  followed  by  elimination  of 
H2NNO2  (Fig.  5).  This  decomposition  path  was  proposed  in  the  literature1  based  on  our 
experimental  products  analysis.  Alternatively,  DNB  could  first  decompose  to  O2N-NHCO- 
NHCO  and  HNNO2,  then  subsequently  undergo  an  inter-molecular  hydrogen  atom  transfer 
between  two  fragments  to  produce  nitramide.  However,  this  path  involves  breaking  a  CN  bond 
with  an  energy  barrier  of  at  least  3.90  eV,  and  thus  is  energetically  unfavorable. 
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Fig.  5:  A  representative  trajectory  of  dissociation  Path  3  at  4000  K,  calculated  at  B3LYP/6- 
31G*.  Plot  shows  two  activated  steps  and  the  variation  of  potential  energy  during  the 
trajectory. 


A  transition-state-theory-based  statistical  model,  Rice-Ramsperager-Kassel-Marcus  (RRKM) 
theory  was  applied  to  examine  the  properties  of  the  reactant  and  transition  states  as  a  function  of 
the  internal  temperature.  RRKM  dissociation  rates  for  initial  decomposition  for  the  various  paths 
are  shown  in  Table  2  for  temperatures  <  1500  K.  RRKM  calculations  predict  that,  at 
temperatures  below  1000  K,  Paths  1  and  2  are  negligible  due  to  the  high  and  tight  transition 
states;  only  formation  of  H2N-NO  would  be  expected.  At  1500  K  (i.e.,  the  maximum  DNB 
denotation  temperature),  decomposition  is  still  dominated  by  Path  3  but  the  contribution  of  Path 
1  becomes  significant.  Note  that  the  calculated  TS3  corresponds  to  the  first  step  of 
decomposition  Path  3,  i.e.,  twist  of  the  O2N-NH  unit  before  H  transfer  (see  Fig.  2  and  Table  1). 
We  were  unable  to  locate  other  TSs  associated  with  the  following  intra-molecular  H  transfer. 
Therefore,  the  RRKM  predicted  reaction  rate  for  Path  3  is  an  upper  limit.  Assuming  that  these 
three  decomposition  paths  account  for  most  DNB  dissociation  at  temperatures  less  than  1500  K, 
the  product  branching  ratios  can  be  estimated  using  the  RRKM  rates.  It  turns  out  that  Path  3 
accounts  for  more  than  80%  of  DNB  decomposition  at  1500  K. 


Table  2.  Calculated  RRKM  dissociation  rates  and  branching  ratios  for  DNB* 


Temperature 

500  K 

1000  K 

1500  K 

k(s-J) 

branching 

k(s‘) 

branching 

k(s‘) 

branching 

Path  1 

0 

0 

9.2  x  108 

0.03 

2.8  x  1010 

0.18 

Path  2 

0 

0 

1.0  x  105 

~0 

5.5  x  108 

~0 

Path  3 

5.1  x  107 

0.01 

2.6  x  1010 

0.07 

1.3  x  10u 

0.82 

The  branching  ratios  were  estimated  when  only  three  major  dissociation  paths  were  included. 
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The  present  work  demonstrates  that  dynamics  simulations,  in  conjunction  with  electronic 
structure  and  RRKM  calculations,  are  able  to  help  extract  decomposition  mechanisms  and 
dynamics  of  complex  energetic  molecules,  product  branching  ratios,  and  how  they  vary  with 
temperature.  The  trajectories  have  been  found  to  reproduce  the  major  decomposition  products  of 
DNB  observed  in  the  pyrolysis  experiment.14  The  trajectories  revealed  many  different 
dissociation  paths.  According  to  RRKM  analysis,  at  least  two  of  these  paths,  i.e.,  elimination  of 
NO2  and  formation  of  nitramide,  are  important  in  the  denotation/hypergolic  experiments.  In 
conclusion,  these  results,  several  of  which  were  difficult  to  anticipate  in  advance,  demonstrate 
the  utility  of  direct  dynamics  trajectory  simulations  as  a  useful  and  necessary  step  for 
investigating  the  decomposition  mechanism  and  dynamics  of  complex  molecular  dissociation, 
and  constructing  kinetics  models.15 

For  the  future,  it  would  be  useful  to  carry  out  further  trajectory  simulations  on  the  subsequent 
dissociation  of  remaining  units  (such  as  O2N-NHCO-NHCO-NH  and  O2N-NHCO-NCO) 
generated  in  the  initial  decomposition  in  order  to  get  a  complete  picture  of  DNB  decomposition. 


2.2  Chemical  Kinetics  Modeling 

Our  FTIR  spectroscopic  measurements1  in  dicyanamide-based  ILs  and  nitric  acid  systems 
showed  the  formation  of  HNCO,  HNO3,  N2O,  NO2,  CO2  and  H2O  as  major  pre-ignition  products. 
Thermodynamic  equilibrium  calculations  on  various  mixtures  of  these  species  were  carried  out 
to  verify  that  such  mixtures  can  sustain  a  flame  under  ambient  experimental  conditions.  To 
further  investigate  which  chemical  pathways  control  the  hypergolic  ignition  process  in  these 
mixtures,  the  thermal  decomposition  of  HNCO,  HNO3  and  N2O  was  written  in  detail  to  explicitly 
account  for  all  of  the  subsequent  reactions  of  the  species  produced  and  their  corresponding  rate 
constants  specified  from  available  thermo-kinetics  databases.  The  detailed  chemical  kinetics 
model  constructed  consisted  of  61  species  and  377  reactions  and  was  implemented  into 
SENKIN.16  The  fate  of  various  homogeneously  premixed  HNCO/CO2/H2O/N2O/HNO3/NO2 
gaseous  mixtures  was  simulated.  The  model  was  able  to  simulate  ignition  at  relatively  low 
temperatures.8  As  shown  in  Fig.  6,  for  a  mixture  at  a  constant  volume  and  at  an  initial 
temperature  of  480  K  and  initial  pressure  of  1  atmosphere,  an  ignition  delay  (ID)  of  19  s  was 
predicted.  Here,  ignition  is  numerically  defined  as  the  time  delay  in  the  system  at  which  point 
the  temperature  profile  has  reached  a  value  of  50%  of  the  maximum.  For  the  composition 
considered  in  Fig.  6,  sensitivity  analysis  showed  that  the  presence  of  NO2  in  the  system  had  a 
significant  effect  on  the  predicted  value  of  the  ignition  delay.  Brute-force  sensitivity  analysis 
suggested  that  the  reaction,  HNO3  — >  OH  +  NO2,  is  an  important  source  for  NO2  in  the  system. 
Typically,  detailed  chemical  kinetics  models  are  validated  against  measured  experimental 
parameters  such  as  ignition  delays,  species  profiles,  and  fundamental  flame  speeds  that  are 
obtained  under  “ideal”  laboratory  setups  with  homogeneous  mixtures  at  known  conditions  of 
temperature  and  pressure.  However,  for  hypergolic  systems,  as  is  the  case  here,  it  is  very 
difficult  to  determine  the  precise  experimental  composition,  pressure  and  temperature  of  the  gas 
phase  in  which  ignition  is  observed  due  to  the  inherent  inhomogeneity  of  the  transiently  formed 
mixtures,  and  therefore  it  is  not  easy  to  validate  the  detailed  chemical  kinetics  model  under 
consideration.  Even  in  the  case  for  premixed  homogeneous  systems,  such  a  validation  effort  in 
many  circumstances  is  not  trivial.  The  present  detailed  chemical  kinetics  model  to  predict 
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ignition  delays  in  ILs  is  the  first  of  its  kind.  Although  the  predictions  are  off  by  2  orders  of 
magnitude,  the  present  effort  represents  a  valuable  initial  step  towards  understanding 
hypergolicity  in  IL/HNO3  systems. 
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Fig.  6:  Computed  temperature  profile  against  time  for  the  mixture 
HNCO/NO2/N2O/CO2/H2O/HNO3:  1/1/2/1/1/1.  The  system  is  considered  adiabatic  and  at 
constant  volume.  P  =  1  atm,  T  =  480  K,  ignition  delay:  19  s. 


Future  refinements  to  the  model  will  be  carried  out  through  an  EOARD  grant.  This  will  include  a 
comprehensive  sensitivity  analysis  both  to  construct  a  reduced-kinetics  model  as  well  as  to 
include  experimentally  determined  new  reaction  paths  and  mechanistic  information  regarding 
NO2  kinetics.  Methods  to  implement  the  reduced-kinetics  model  into  a  computational  fluid  code, 
such  as  FLUENT,  will  be  explored.  Techniques  to  validate  the  resulting  model  against  realistic 
experimental  parameters  that  describe  transient  composition,  pressure,  and  temperature  profdes 
during  hypergolic  ignition  of  IL/HNO3  mixtures  will  be  assessed.  The  ultimate  capability  we 
seek  in  this  modeling  and  simulation  effort  is  the  ability  to  accurately  predict  ID  times  in  known 
IL/oxidizer  combinations  as  well  as  to  be  able  to  numerically  identify  hitherto  unknown  hypergol 
combinations. 


3.0  EXPERIMENTAL  EFFORTS 
3.1  FTIR  Spectroscopic  Investigations 

To  understand  the  initiation  processes  that  occur  in  the  condensed  phase  during  the  reaction  of  an 
IL  with  various  oxidizers,  transmission  FTIR  spectroscopic  investigations  were  performed  to 
elucidate  how  the  oxidizer  initially  ionizes  upon  contact  with  the  IL.  Subsequent  reactions  of 
these  ionized  species  with  the  ionic  components  of  the  IL  were  followed  by  recording  the 
temporal  behavior  of  the  IR  absorptions  of  the  IL  in  the  condensed  phase. 
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Experiments  were  performed  by  preparing  a  thin- film  of  the  IL  on  a  Ge  disc  and  exposing  it  to 
an  excess  of  the  oxidizing  gas  under  anhydrous  conditions.  Simple  ionization  of  both  NO2  (<-> 
V2N2O4)  and  HNO3,  respectively,  to  NO+  +  NO3  and  H+  +  NO3  was  observed  to  occur  at  the 
ionic  surface  as  seen  in  Fig.  7.  Transient  formation  of  NO+  was  inferred  from  FTIR  observation 
of  the  N2O  product  when  the  ionic  liquid,  R'Nj"  was  treated  with  N2O4,  where  the  reaction  N3'  + 
NO+  leads  to  N20  +  N2  formation. 

The  difference  in  the  heterogeneous  chemistry  of  NO2  and  HNO3  towards  a  thin- film  of  1-butyl- 
3-methyl-imidazolium  dicyanamide  (BMIM'N(CN)2~)  is  shown,  respectively,  in  Figs.  8  and  9. 
The  spectrum  of  unreacted  BMIM+N(CN)2‘  (blue  trace,  Figure  8)  is  shown  with  characteristic 
N(CN)2-  fundamental  and  combination  bands  identified.  These  bands  disappear  upon  reaction  of 
BMIM+N(CN)2_  with  gaseous  NO2  (red  trace,  Figure  8),  providing  direct  spectroscopic 
confirmation  that  the  anion  initiates  the  reaction  with  the  oxidizing  species.  All  of  the  cationic 
bands  in  the  blue  trace  are  also  present  in  the  red  trace,  and  in  addition,  the  appearance  of  new 
product  bands  is  noted.  This  suggests  that  the  IF-cation  merely  acts  as  a  spectator  during 
initiation,  until  after  ignition  has  occurred,  whereupon  it  would  then  participate  in  further  high 
temperature  oxidative  processes  during  steady-state  combustion  of  the  IL. 
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Fig.  7:  NO3'  absorption  bands  in  the  ionization  of  NO2  (red)  and  HNO3  (black)  at  an  IL 
surface. 


The  product  spectrum  of  BMIM+N(CN)2_  after  exposure  to  gaseous  HNO3  (purple  trace  in  Fig.  9) 
also  indicates  that  the  anion  is  initially  involved  in  the  reaction.  However,  different  product 
bands  are  noted  in  this  trace  compared  to  that  in  the  red  trace  of  Fig.  8.  Comparison  to  the 
experimentally  recorded  biuret  spectrum  suggests  that  amide-III  like  features  (peaks  19  and  20) 
are  also  present  in  this  product  spectrum  (peaks  32,  35,  and  36).  Characteristic  features  of  N-H 
stretching  in  the  purple  trace  were  also  noted. 
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Fig.  8:  Spectrum  of  neat  thin-film  of 
BMIM+N(CN)2'  (blue  trace)  and  after 
reaction  with  N02  (red  trace). 


wavenumber 

Fig.  9:  Product  spectrum  of  thin-film 
of  BMIM+N(CN)2‘  after  reaction  with 
gaseous  HNO3  (purple  trace)  and  of 
neat  biuret  (green  trace). 


To  confirm  the  possibility  that  the  1,5-dinitrobiuret  (DNB)  intermediate  is  being  formed  and 
trapped  in  the  thin-film,  ab  initio  computations  of  the  DNB  infrared  spectrum  as  well  as  of  other 
possible  intermediates  in  BMIM+N(CN)27N02  and  BMIM+N(CN)27HN03  are  currently 
underway.17  Our  preliminary  analysis  of  computed  spectra  at  the  B3LYP/6-31+G(d,p)  level  of 
theory  are  consistent  with  DNB  formation. 

Our  future  goal  in  the  experimental  FTIR  spectroscopy  is  to  significantly  expand  our 
understanding  of  the  heterogeneous  phase  chemistry  through  detailed  kinetics  measurements  of 
product  and  reactant  temporal  profiles  at  ionic  liquid-oxidizer  gas  interfaces.  We  plan  to  do  this 
by  employing  an  attenuated  total  reflection  (ATR)  cell  in  which  there  will  be  temperature  control 
of  the  ionic  liquid  thin-film  as  well  as  pressure  and  concentration  control  of  the  oxidizer.  All 
three  oxidizers,  HNO3,  N2C>4  and  H202  will  be  investigated.  This  will  allow  reaction  orders  to  be 
calculated  as  well  as  Arrhenius  parameters  for  the  products  formed  in  each  case.  The  ATR  cell 
will  provide  a  capability  to  elucidate  the  effect  of  impurities  on  initiation  processes  in  the  IL. 

This  information  is  of  particular  importance  in  ignition  model  development;  since  often  the 
purity  of  ionic  liquid  fuel  can  vary  considerably,  which  can  significantly  affect  its  ability  to 
hypergolically  ignite.  For  instance,  the  influence  of  water  impurity  will  be  investigated  by 
introducing  a  known  amount  of  water  both  directly  into  the  ionic  liquid  as  well  as  an  added 
vapor  component  in  the  oxidizer  gas. 

3.2  Ignition  Delay  Measurements 

The  search  for  new  igniting  room  temperature  IL/oxidizer  combinations  continued  with  the 
discovery  of  substituted  imidazolium  thiocyanate-based  ILs  and  HNO3  hypergols  using  our  fast- 
videography  ignition  delay  (ID)  apparatus.1  Future  studies  of  the  reaction  chemistry  of  the  SCN' 
anion  with  HNO3  would  be  of  interest  to  contrast  the  proposed  initiation  process  in  N(CN)2" 
/HNO3  system. 
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The  role  of  complex  ionization  of  the  HNO3  when  mixed  with  various  ionic  liquids  is  also  being 
investigated.  The  extent  of  nitronium  ion  (N02+)  formation  in  the  process:  2H+  +  2NO3' 

<=>  2HNO3  <=>  N02+  +  H20  +  NO3'  can  be  controlled  by  the  addition  of  excess  H2SO4  to  the  nitric 
acid:  2H2SO4  +  HNO3  <=>  N02+  +  H30+  +  2HSCV.  The  effect  of  [N02+]  on  both  the  ignition 
delay  as  well  as  the  composition  of  the  pre-ignition  products  in  BMIM+N(CN)2"/H2S04/HN03 
mixtures  is  currently  being  studied  using  our  rapid-scan  FTIR  machine. 

3.3  Vacuum  Ultraviolet  (VUV)  Spectroscopic  Investigations 

The  above  work  was  further  complemented  by  carrying  out  more  detailed  investigations  of  room 
temperature  IL  vaporization  and  pyrolysis  processes.  A  thorough  understanding  of  the  relative 
role  of  vaporization,  pyrolysis  and  oxidation  when  ILs  hypergolically  react  with  HNO3  is 
important  for  constructing  high  fidelity  kinetics  models  to  predict  ID  times.  Vaporization  studies 
were  performed  at  the  Chemical  Dynamics  Beamline  at  the  Advanced  Light  Source  (ALS) 
facility  of  the  Lawrence  Berkeley  National  Laboratory  in  collaboration  with  Professor  Steve 
Leone.  This  work  was  carried  out  under  ALS  beam-access  Grant  #  02756.  The  Chemical 
Dynamics  Beamline  at  the  ALS  is  capable  of  producing  high-intensity  continuous  vacuum 
ultraviolet  (VUV)  light  in  the  range  of  7.4-15  eV,  which,  when  coupled  to  a  time  of  flight  mass 
spectrometer,  is  appropriate  for  determining  photoionization  efficiency  (PIE)  curves.  PIE  data 
allows  for  the  determination  of  photoion  appearance  energies  and  ionization  potentials  which 
provide  a  powerful  means  to  identify  the  gas-phase  species  present  when  ILs  are  heated  in 
vacuum. 

Experiments  were  carried  out  to  identify  the  vapors  evolved  when  ILs  are  heated.  Since  ILs 
have  extremely  low  vapor  pressures,  identification  of  the  vapors  evolved  has  been  difficult  to 
perform  experimentally.  Our  research  indicates  there  exist  three  regimes  important  in  ionic 
liquid  vaporization,  depending  on  the  IL  heats  of  vaporization  (AHvap)  relative  to  their  lowest 
activation  barrier  process  to  thermal  decomposition  (Ea):  1)  AHvap  «  Ea;  2)  AHvap  ~  Ea;  and  3) 
AHvap  »  Ea.  In  Case  1,  the  vapors  are  solely  from  vaporization  of  the  ionic  liquid,  typically  as 
intact  cation-anion  ion  pairs.2,3  In  Case  3,  since  the  heat  of  vaporization  is  much  higher  than  the 
barrier  to  thermal  decomposition,  the  vapors  will  consist  exclusively  of  the  volatile  IL  thermal 
decomposition  products.  In  Case  2,  there  will  be  a  competition  between  vaporization  of  the  IL 
and  thermal  decomposition  of  the  IL,  and  the  vapor  will  consist  of  a  combination  of  vaporized  IL 
and  volatile  IL  thermal  decomposition  products.  Examples  of  each  case  are  as  follows:  1) 
EMIM'NTf2";  2)  EMIM+Br';  and  3)  EMIM+SCN~.  The  nucleophilicity  of  the  anion  could  also 

1  a 

play  an  important  role  in  the  thermal  stability  of  these  ILs. 

Our  ALS  studies  have  demonstrated:  1)  the  influence  of  molecular  structure  on  the  photoion 
appearance  energies  of  isolated  ion  pairs;  2)  heats  of  vaporization  of  ILs  can  be  determined  via 
VUV  time-of-flight-mass-spectrometry  (TOFMS);  3)  the  nature  of  thermal  decomposition 
pathways  of  several  ionic  liquids;  and  4)  the  occurrence  of  new  dissociative  photoionization 
processes  for  some  gas  phase  ionic  liquids.  Also,  studies  with  aprotic  ILs  containing  strongly 
nucleophilic  anions  (N3',  Br,  N(CN)2",  SCN')  indicate  a  variety  of  reactions  occurring  when 
heated:  polymerization,  proton  or  alkyl  transfer  to  the  anion,  and  addition  of  the  anion  to  the 
cation  followed  by  molecular  elimination. 
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The  extremely  low  vapor  pressures  of  room  temperature  ILs  present  a  challenge  to  experimental 
design  for  studying  their  reactivity,  as  they  are  not  conducive  to  mass  spectrometric  study. 
However,  by  generating  IL  aerosols  we  can  circumvent  this  issue.  The  process  by  which 
gaseous  IL  species,  in  the  form  of  aerosolized  particles,  oxidize  in  HNO3  is  currently  the  subject 
of  investigations  under  a  new  ALS  beam-access  Grant  #  03341.  Recently,  generation  of  ionic 
liquid  aerosols  (see  Fig.  10)  has  been  achieved  by  standard  aerosol  techniques.19  Optimization 
studies  of  the  highly  thermally  stable  ionic  liquid,  EMIM+NTf2',  indicates  that  0.5  g/1  of  IL  in 
water  is  a  good  choice  for  producing  -150  nm  IL  aerosols  which  are  expected  to  be  essentially 
solvent  free.  Less  thermally  stable  ILs  relevant  to  Air  Force  interests  have  been  shown  to  form 
aerosols  as  well.  Preliminary  work  is  underway  at  the  Air  Force  Research  Laboratory  to 
determine  safe  operating  conditions  for  reacting  BMIM+N(CN)2'  aerosols  with  various  oxidizing 
vapors. 


Ionic  Liquid  Aerosol 

Emim"Tf2N'  (in  water) 

CHi 

O  O 

£  S  F^C-S-N-S-CFj 
N  6  6 


Cone.  Peak  diameter 


Mean  diameter 


0.2g/L 

25.0  nm 

31.6  nm 

0.5g/L 

34.6  nm 

50.5  nm 

1.0g/L 

41.4  nm 

50.8  nm 

Fig.  10:  Size  distributions  for  ionic  liquid  aerosols  generated  from  varying  concentrations 
in  water. 


These  experimental  conditions  will  be  reproduced  at  the  aerosol  end-station  available  at  ALS 
Beamline  9.0.2.  Initial  experiments  will  characterize  the  IL  aerosols,  including  aerosol  size 
distributions,  remaining  solvent  content,  and  evaluate  the  photoionization  mass  spectra  in  the 
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absence  of  oxidizer  to  establish  a  baseline.  Then,  the  dicyanamide-based  room  temperature  IL 
aerosols  will  be  reacted  under  dilute  HNO3  conditions  and  the  PIE  of  the  gaseous  products  will 
be  measured.  Finally,  analysis  of  the  reacted  aerosol  particles,  where  the  particles  are  impinged 
on  a  hot  surface  to  vaporize  the  products,  will  allow  for  the  identification  of  condensed  phase 
products.  Comparison  of  oxidation  products  with  those  in  the  proposed  pre-ignition  mechanism1 
will  be  carried  out.  Expected  products  include  HNCO  (mass  43,  IP=1 1.6  eV),  N2O  (mass  44, 
IP=12.9  eV),  C02  (mass  44,  IP=13.8  eY),  NH3  (mass  17,  IP=10.07  eY),  and  H2NN02  (mass  62). 
By  fixing  the  photoionization  energy  and  varying  other  experimental  parameters  such  as 
temperature,  oxidizer  concentration,  aerosol  particle  diameter,  and  mixing  time  in  the  flow  tube, 
kinetics  measurements  will  be  performed  to  determine  the  reactive  uptake  coefficients  and 
activation  energies  involved.  The  investigations  will  be  extended  to  alternative  oxidizers  such  as 
NO2  and  H2O2. 


4.0  FUTURE  OUTLOOK 

Upon  successfully  carrying  out  the  oxidation  reactions  on  the  dicyanamide-based  ionic  liquids, 
we  will  systematically  extend  our  studies  to  new  classes  of  high  performance  hypergolic  room 
temperature  ILs  where  the  pre-ignition  mechanisms  are  not  known  or  have  not  been  investigated 
experimentally.  Elucidation  of  the  oxidative  mechanisms  in  these  new  compounds  and 
comparison  to  the  dicyanamide-based  room  temperature  ILs  will  provide  valuable  information 
regarding  their  dramatically  different  ignition  characteristics  that  we  have  recently  observed  at 
the  Air  Force  Research  Laboratory  (AFRL).  The  development  of  fast  chemical  (hypergolic) 
ignition  methods  for  these  room  temperature  ILs  is  of  high  importance  to  the  Air  Force  for 
designing  new  energetic  bipropellants  for  space  propulsion  applications.  A  detailed 
understanding  of  the  oxidation  and  pre-ignition  chemistry  that  will  be  possible  at  the  AFRL 
Propellants  Branch  and  ALS  facilities,  along  with  the  previously  studied  competing  thermal 
decomposition  processes  that  are  inevitable  as  the  temperature  of  the  room  temperature  IL  is 
raised  during  hypergolic  combustion,  will,  for  the  first  time,  allow  us  to  construct  chemical 
kinetics  models  to  predict  the  ignition  behavior  of  these  fuels.  These  proposed  laboratory 
aerosol-oxidizer  studies  are  one  step  closer  to  real-world  conditions  within  a  propulsion  system, 
where  IL  fuels  are  better  represented  by  sprays  or  liquid  droplets  mixing  with  vaporized  oxidizer 
under  high  temperatures  and  pressures. 


5.0  CONCLUSIONS 

New  hypergolic  IL  fuels  continue  to  be  discovered  by  the  trial-and-error  method  of  first 
synthesis  and  then  test  approach.20'22  Our  desire  to  understand  the  fundamental  chemistry  of 
initiation  and  ignition  in  IL  hypergols  provides  an  opportunity  to  construct  numerical  models  for 
designing  fuels  which  have  the  required  reactivity  using  “first  principle”  theoretical  tools  prior  to 
synthesis.  We  have  seen  preliminary  success  of  our  efforts  in  a  chemical  kinetics  model,8 
quantitative-structure-property-relationship  correlation  model,10  and  MD  simulations.9,15 
Continued  investigations  to  understand  the  interplay  between  physical  phenomena,  i.e.,  mixing 
and  interfacial  transport,  and  chemical  kinetics,  i.e.,  key  initiation  reactions,  pre-ignition 
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intermediates  and  the  ignition  event,  will  greatly  facilitate  the  efficient  discovery  of  new 
hypergols  of  importance  to  the  IHPRPT  Program. 
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A  class  of  room- temperature  ionic  liquids  (RTILs)  that  exhibit  hyper  go  lie  activity  toward  fuming  nitric  acid 
is  reported.  Fast  ignition  of  dicyanamide  ionic  liquids  when  mixed  with  nitric  acid  is  contrasted  with  the 
reactivity  of  the  ionic  liquid  a/iti&s,  which  show  high  reactivity  with  nitric  acid,  hut  do  not  ignite.  The  reactivity 
Of  other  potential  sail  fuels  is  assessed  here.  Rapid-scan,  Fourier  transform  infrared  (FTIR)  spectroscopy  of 
the  preignition  phase  indicates  the  evolution  of  N20  from  both  the  dicyanamide  and  azide  RTELs,  Evidence 
for  the  evolution  of  CCh  and  isocyanic  add  (HNCO)  with  similar  temporal  behavior  to  NT)  from  reaction  of 
the  dicyanamide  ionic  liquids  with  nitric  acid  is  presented.  Evolution  of  HN:*  is  detected  from  the  azides.  No 
evolution  of  i  ICN  from  the  dicyanamide  reactions  was  detected.  From  the  FTIR  observations,  hi urd  reaction 
tests,  and  initial  ab  initio  calculations,  a  mechanism  is  proposed  for  the  formation  of  N?Of  CO?,  and  HNCO 
from  the  dicyanamide  reactions  during  preignition. 


Introduction 

There  continues  to  be  a  growing  effort  to  synthesize  task- 
specific  ionic  liquids  (M.s).  Ionic  liquids  arc  sails  wilh  melting 
points  below  100  °C.  Lately,  the  design  and  choice  of  many 
ionic  liquids  is  focused  on  physical  properties  such  as  miscibil¬ 
ity,  conductivity,  viscosity,  solubility,  and  melting  points.1  An 
interesting  and  important  property  of  ionic  liquids  is  Thai  they 
have  very  low  vapor  pressures,  with  high  heats  of  vaporization 
(approaching  ^60  kcal  mol  1)r  Tire  details  of  how  the  chemical 
structure  of  the  ionic  liquid  nf  fee  I  s  these  various  physical 
characteristics  are  now  beginning  Lo  emerge/1  ‘ 6  Less  understood, 
but  now  of  great  interest,  arc  the  associated  chemical  properties 
that  manifest  in  the  condensed  phase  such  as  relative  basicity 
(or  acidity),  oxidative  (or  reductive)  capacity,  then  rial  stability, 
electrochemical  reactivity,  and  catalytic  or  combustion  efficiency 
of  the  component  ions.7  9 

Recent  interest  in  room  temperature  ionic  liquids  has  grown 
immensely,  as  their  synthetic  routes  have  been  optimized.  The 
number  of  publications  involving  ionic  liquid  chemistry  has 
grown  exponentially  in  the  past  decade.1-10-11  '[be  number  of 
possible  different  cation— anion  combinations  to  form  ILs  is 
enormous,  and  so  it  is  important  to  understand  trends  hi 
reactivity  and  to  have  theoretical  tools  to  screen  for  the  best  IL 
candidates  for  a  given  use.  Typical  cations  for  ionic  liquids 
considered  here  include  asymmetrically  substituted  imidaze  limns 
and  triazoliums,  since  symmetric  cations  tend  to  produce  solid 
salts  at  room  rem.perarure,6  However,  for  a  given  cationic  system, 
the  choice  of  the  anion  is  also  important  in  tailoring  die  liquidus 
range  of  die  salt, 12  Similarly,  the  chemical  and  thermal  stabilities 
of  the  resulting  !L  can  depend  both  on  Ihe  physical  (structural) 
and  chemical  properties  of  Ihe  individual  component  ionsA15,3'1 
For  instance,  in  the  case  of  hnidazolium  ionic  liquids,  studies 
have  shown  that  the  overall  thermal  stability  of  the  IL  correlates 
with  the  anion's  size,™  its  micleoph  ill  city  (or  its  Lewis 
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basicity),13-16’17  and  its  liydrophilicity,1*  etc.  On  the  other  hand, 
for  a  given  anionic  system,  cationic  si/e  has  little  effect  on  Ihe 
IL  stability,™1*  w^hile  increased  substituiion  of  die  imidazolium 
ring  at  die  C(2) -site  by  straight- alkyl  chains  increases  the 
stability  but  secondary  or  tertiary  alkyl  groups  at  the  N-site 
decreases  Ihe  stability.15 

Recently,  several  groups  have  shown  that  some  classes  of 
ionic  liquids  can  be  distilled  in.  a  vacuum,20  These  ILs  tend  to 
have  aprotic  cations  and  non-uncloophilic  anions,  and  studies 
have  shown  that  the  gas -phase  species  exist  (in  high  vacuum) 
as  neutral  ion  pairs  and  do  not  tend  to  vaporize  as  gas -phase 
clusters.21”25  For  pi  otic  ionic  liquids,  however,  die  formation 
of  neutral  molecules  in  the  gas  phase  has  also  been  reported,22 
An  important  Ihmnucha ideal  properly  which  is  difficult  to 
obtain  accurately  for  ionic  liquids  is  die  condensed  phase  heal 
of  formation.  Due  to  the  ionic  nature  and  liquid  phase  of  ILs, 
theoretical  calculations  of  the  condensed  phase  heats  of  forma¬ 
tion  are  d  if  tic  ul  r  and  have  large  uncertainties  associated  with 
them.  Recent  experimental  measurements  of  heats  of  combustion 
and  heats  of  vaporization  of  ILs  have  begun  to  accurately 
determine  the  condensed  phase  heat  of  formation  for  a  few  ionic 
liquids.2^ 

Here,  we  present  die  results  of  our  efforts  lo  tailor  ionic 
liquids  for  utilization  in  bipropellant  rocket  engines.  We  desire 
not  only  to  improve  the  performance  (he.,  density  specific 
impulse)  over  ihe  current,  stule-oJ-the-arl,  monomethyi  hydra¬ 
zine/nitrogen  tetroxide  (MMH/NTG)  system,  but  also  be  able 
to  increase  die  ambient  reactivity  of  the  component  ions  toward 
common  oxidizers,  such  that  hypergolic  ignition  (spontaneous 
ignition  upon  mixing)  is  attained,  comparable  to  MMII/NTO 
ignition  times.  To  calculate  performance  characteristics  of 
propellants,  the  condensed  phase  heat  of  formation  value  is 
critical.2^27  To  successfully  implement  such  ionic  liquids  into 
viable  working  systems,  their  physical  properties  such  as 
viscosity  and  the  liquid  temperature  range  must  also  be  tailored 
within  required  operational  ranges. 1*  Furthermore,  the  inherent 
low  vapor  pressure  of  ionic  liquids  should  significantly  reduce 
the  ambient  toxicity  and  environmental  emissions  and  therefore 
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Figure  |t  Schematic  of  the  FTIR  apparatus.  The  inset  illustration  in  the  lower  left  comer  shows  the  relative  positions  of  the  HeNe  laser  beam,  the 
IR  beam,  and  the  falling  oxidizer  droplet. 


the  associated  cost  of  handling  and  fielding  of  these  new 
replacement  fuels.  The  extremely  low  vapor  pressure  of  ionic 
liquids  at  ambient  temperature  and  pressure  is  a  major  benefit 
of  hypergolic  ionic  liquid  fuels  versus  MMII,  which  is  highly 
volatile.  The  high  vapor  toxicity  of  MMII  increases  the  safely 
precautions  required,  and  the  logistics  of  handling  i  t  safely  while 
loading  ihe  space- vehicle  can  he  very  expensive. 

Understanding  of  chemical  pathways  and  reaction  mecha¬ 
nisms  when  fuels/]  Ls  react  with  strong  oxidizers  is  critical  to 
designing  new  fuels.  The  current  MMH/NTO  hypergolic  system 
has  been  used  successfully  for  decades,2*-29  and  much  theoretical 
work  on  the  MMH/NTO  system  has  been  performed,30-35 
However,  little  experimental  work  exists36-40  to  confirm  these 
results.  Therefore,  replacing  the  volatile,  highly  toxic,  and 
carcinogenic  MMII  with  a  safer  fuel  is  not  a  straightforward 
task.  One  of  the  challenges  is  I  hat  the  mechanism  involves 
highly  reactive  transient  species  which  are  difficult  to  probe 
experimentally.41,42  The  identification  of  common  reactive 
species  in  different  hypergolic  systems  will  play  a  significant 
role  in  future  fuel  design.  There  are  two  different  regimes 
involved  in  hypergolic  systems:  first,  the  low-temperature 
preignilion  phase,  where  gaseous  reactive  species  are  produced 
and  the  temperature  of  the  system  increases  as  thermolytic  and 
oxidative  processes  in  the  condensed  phase  begin  to  decompose 
both  the  fuel  and  the  oxidizer.  Once  the  temperature  and 
concentrations  of  key  species  reach  a  critical  level,  ignition 
occurs  in  the  gas  phase,38-43  the  temperature  increases  dramati¬ 
cally,  and  combustion  occurs.  To  construct  accurate  numerical 
combustion  models  for  IL  hypergols,  detailed  chemical  mecha¬ 
nistic  information  is  needed  on  how  fuel  vaporization/pyrolysis 
processes  compete  with  and  augment  concurrent  and  subsequent 
oxidative  processes  of  the  reactive  fragments  that  are  being 
produced. 

Initial  attempts  by  the  Air  Force  Research  Laboratory  to 
design  ILs  with  highly  strained  or  unsaiuraled  functional  groups 
were  made  to  introduce  '"trigger  groups"  onto  the  cation  in  the 
hope  of  inducing  hypergolicity.  These  cations  may  he  paired 
with  oxygen-containing  anions  (NO^-)  for  oxygen  balance  or 
high  heals  of  formation  (N3  )  to  promote  hypergolic  ignition 
when  treated  with  a  suitable  oxidizer.44  The  discovery  that  the 
anion  can  play  a  critical  role  in  inducing  hypergolicity  Is 
described  here.  The  fact  that  a  common  anion  can  induce 
hypergolicity  when  paired  with  numerous  different  cations  to 
form  ionic  liquids  allows  for  the  tailoring  of  the  cation  for 


various  desirable  properties  such  as  high  heat  of  formation,  low 
viscosity,  a  wide  liquid  temperature  range,  environmental  safely, 
and  thermal  and  shock  stability.  This  paper  describes  the 
investigation  of  the  preignition  phase  of  hypergolic  ignition  of 
ionic  liquids,  a  critical  anion  needed  to  induce  hypergolic 
reactivity,  and  the  influence  of  cation  structure  on  the  reactivity. 
This  investigation  focuses  on  white  fuming  nitric  acid  (WFNA) 
as  the  oxidizer, 

K  x  peri  mental  T celt  m  q  u  e 

Fouricr-transform  infrared  (FTIR)  absorption  spectroscopy 
was  performed  using  a  rapid-scan  spectrometer.  Typical  spectra 
were  acquired  in  the  range  590—3850  cm  1  with  a  spectral 
resolution  of  2  cm-1  and  approximately  every  30  ms.  The 
interferograms  were  recorded  in  split  double- sided,  split 
forward— backward  mode.  The  reaction  was  carried  out  within 
a  closed  stainless  steel  chamber  under  constant  nitrogen  purge 
and  therefore  was  carried  out  in  an  oxygen- free  environment. 
Figure  1  shows  the  schematic  of  our  experimental  setup.  The 
infrared  beam  was  focused  into  the  chamber  above  the  reaction 
cell  by  means  of  gold  minors  and  KRS-5  windows.  A  drop  of 
oxidizer  was  introduced  from  a  gaslight  syringe  through  a 
septum  at  the  lop  of  the  chamber.  The  drop  fell  into  a  small 
cuvette  containing  a  small  amount  (0,1  —  1 .0  mL)  of  ionic  liquid 
fuel.  The  spectrometer  was  triggered  via  a  IleNc  laser  and 
photodiode  positioned  above  the  reaction  cell.  As  the  drop  of 
oxidizer  passes  through  the  IIcNc  beam,  the  drop  in  signal  from 
the  photodiode  triggers  the  spectrometer.  Data  was  collected 
before  the  ignition  flash  occurs  and  is  hereafter  referred  to  as 
the  preignilion  phase  data. 

To  observe  and  measure  the  actual  ignition  delay  (ID),  a  small 
fuel  sample  (10— 50  h«L>  was  instead  dropped  into  the  cuvette 
containing  an  excess  amount  (1 ,0  mL)  of  the  oxidizer,  A  high¬ 
speed,  complementary  metal-oxide  semiconductor  (CMOS) 
video  imager  was  used  to  record  the  time  duration  from  the 
moment  the  fuel  comes  in  contact  with  the  liquid  surface  of 
the  oxidizer  and  the  first  sign  when  a  visible  flame  is  formed, 
which  invariably  was  always  in  the  gas  phase  above  the  liquid 
surface.  In  the  present  work,  the  CMOS  video  settings  of  1000 
frame s/see  and  168  x  512  pixel  resolution  were  chosen.  Figure 
2  depicts  typical  observations  for  MMII/NTO,  MMH/WFNA, 
and  l-propargyl-3-methyl-imidazolium  dicyanamide/WFNA  hy¬ 
pergolic  ignitions  for  which  the  ID  times  were  determined  to 
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a  b  c 


Figure  2,  Hypergolic  ignition  of  (a)  MMH/NTO,  ID  =  1  ms;  (b)  MMH/WFNAt  ID  =  9  m$;  and  (c)  l-pmpargyl-3-methyl-imida2olium  dieyanamide/ 
WFNA.  ID  =  15  ms. 

TABLE  t;  Theoretical  Reaction  Enthalpies  Calculated  at  (he  Hartrec— Fock  and  Density  Functional  Theory  Levels 

reaction 

Af/oK  kcal/mol  (zpe  corrected) 

method 

NOr  +  H4  —  IlNOj 

—318.5 

B3-LYP  6-3  lG(+)(d,p) 

N(CN)2“  +  IT  —  HN{CM)i 

—302*8 

B3-LYP  6-3lG(+Xd.p) 

N(CN)a  +  H4  —  (NC)N(CNH) 

-310*4 

B3-LYP  6-31G(+Xd,p) 

11N(CN):  +  ir  —  lJ.2N(CN)a+ 

-130.9 

B3-LYP  6-31G(+Xd,p) 

(NC)N(CNH)  +  H1  —  (NC)NH(CNH)4 

-166.2 

B3-LYP  6-3lG(+Xd,p) 

N(NO;):  +  IT  —  l-IN(NO»fc 

-303.8 

RHF  6-31G(+)(d,p) 

N(NO:)r  +  IT  —  1IN(N02)2 

-303.6 

B3-LYP  6-31G(+Xd,p) 

N(N03h  +  H~  —  (02N)N(NO:H) 

-297*0 

RHF6-3lGC+)(d,p) 

HN(N02>2  +  H4  —  H:N(N02)24 

-169,9 

B3-LYP  6-31G{+Xd,p) 

(OaN)N(N0iH)  +  IT  —  H20  +  NaO  +  NO/ 

-21X4 

RHF  6-3lG(+)(d,p) 

HN(NO ih  T  H+  —  H2O  +  N>G  +  NQ24 

-205*5 

RHF  6-31G(+)(d.p) 

(NC)N(N02>-  +  K*  —  (NC)NII(NOa) 

-307,8 

RHF  6-31G(+)<d.p) 

(NC)N(NOi)-  +  iT  —  (HNQN(N02) 

-297*9 

RHF  6AlG(+)(d,p) 

(NC)N(N02)  +  H  *  (NC)N{N02H) 

-321*6 

RHF  6-31G(+)(d.p) 

(NC)N(NOjH)  +  IT  —  (NC)N(N02Ha)+  —  H2O  +  NCNNO4 

-122,6 

RHF  6-3lG(+Kd.p) 

(NC)N(NOzH)  +  H+  —  HaO  +  NaO  +  CN+ 

28*2 

RHF  6-3lG(+)(d,p) 

(NC)NH(NG2)  +  H4  —  HaO  +  NaO  +  CN' 

14,4 

RHF  6-31G(+)(d.p) 

(HNC)N(NOi)  +  H4  —  II3O  +  NaO  +  CN4 

4.5 

RHF  6-31G(+)(d,p) 

(NC)N(NOaH)  +  H+  —  (NC)NH(NOaHj4 

-14B.I 

RHF’  6-3 lG(+)(d,p) 

(NC)N(NOiH)  +  IT  —  (HNG)N(NO:H)+ 

-166,2 

RHF  6-31G(+)(d.p) 

(HNC)N(NOa)  +  H+  -  (HNONIKNO/4 

-18X5 

RHF  6-31G(+)(d.p) 

(NQNH(NOi)  +  IT  —  (HNC)NH(N02)+ 

-172,6 

RHF  6  31G(+)(d.p) 

(HNC)N(NOa)  +  H4  —  (llNQN(NOzH)+ 

-189.9 

RHF6-310(+)Cd,p) 

NC-N-CNH  +  N03  —  NC-N-C(N03}NH 

-5,4 

B3-LYP  6-3lG(+){d.p) 

NC-N-C(NOONH  —  NC-N-C(=0)NfN0:hlT 

-34*5 

B3-LYP  6-31G(+)(d,p) 

NC— N-C(— 0)N{N0;jl  1  +  H4  —  HNC-Nt:(=O)N(N02)H 

-305*2 

B3-LYP  6-3  IG(+Xd,p) 

HNC-N-C(=0)NtN02^  +  H4  —  HNC-NH-C{=O)N(N02)H4 

-178*0 

B3-LYP  6-31Gf+Xd.p) 

be  9  and  15  ins,  respectively.  For  the  interested  reader  we 
provide  Supporting  Information  as  a  mpg-tile  that  shows 
hypergolic  ignition. 

On  the  basis  of  a  common  understanding  of  how  certain 
energetic  materials  are  known  to  decompose  to  stable  products, 
and  the  similarity  in  the  products  we  observed  here,  a  suitable 


method  to  identify  the  nature  of  at  least  one  of  the  possible 
intermediates  was  carried  out  as  follows.  The  reaction  of  ILs 
with  nitric  acid  can  be  slowed  down  using  34  wt  %  HNO3 
solutions  instead  of  WIT4A,  This  allows  for  the  possibility  to 
trap  reaction  intermediates  by  quenching  and  then  suitabiy 
identifying  their  nature,  f  [ere.  we  performed  the  biuret  test  using 
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the  method  of  Lowry45  io  identify  the  presence  of  peptide-like 
functionality  in  the  reaction  intermediate  using  Cu2~  solutions. 
Confirmation  of  the  presence  of  such  a  functionality  together 
with  the  observed  and  theoretical  evaluations  of  all  the 
systematically  substituted  anions  studied  here  would  provide 
strong  evidence  tor  the  proposed  mechanism  (see  Discussion). 

Materials.  The  following  chemicals  were  synthesized  and 
characterized:  l-butyl-3-methyl-imidazoIium  azide  (>97%), 
l-aJlyI-3-mcthyl-imidazolium  azide  (>97%),  l-butyl-3-methyl- 
imidazolium  dieyanamide  (>97%),  l-allyl-3-meihyl-imidazo- 
liurn  dieyanamide  (>97%),  1  -propargyl-3-methyl-i  midazolium 
dieyanamide  (>97%),  l-methyl-4-amino-L2T4T-Lriazolium  di- 
cyanamide  (>97%),  1 -( 3-b u teny l>- 3-methyMmjdazolium  dicy- 
anamide  (>97%),  l-(2-pentynyl)-3-methyl-injidazotium  dicy- 
anamide(>97%),  1  -methyM-amino- 1 ,2,4,-triazr^iumniuxx.yanamidc 
(>97%),  LalJyl-4-amintv l  ,2,4-tnazolium  nitix Ky anamide  ( > 97%), 
and  silver  oitrocy  an  amide  (>97%).  The  characterization  data 
of  the  a  hove  azides  is  reported  elsewhere.44  The  purities  of  the 
above  compounds  were  determined  by  elemental  analysis. 
Potassium  dinilramide  (>97%)  was  kindly  supplied  by  SRI.  The 
following  chemicals  were  purchased  from  Merck:  1-butyl- 3- 
met  hyl-i  mida  zo  I  i  u  m  di  c  y  an  a  mi  de  ( >  99% ),  I  -bmy  I  - 1  -  methy  l- 
pyrolidinium  dieyanamide  (>99%),  and  n-butyl-3-meihyL 
pyridiniu  m  dieyanamide  ( >  99%).  1  -ethyl -3- met  by  l-imidazohuin 
dieyanamide  (>98%)  was  purchased  from  Ruka.  Sodium  azide 
(>99.5%),  sodium  dieyanamide  (>96%),  and  silver  cyanate 
(>99%)  were  purchased  from  Sigma  Aldrich,  and  while  fuming 
nitric  acid  (>99.5%)  was  purchased  from  Fluka.  All  chemicals 
purchased  were  used  withoul  further  purification. 

Ah  Initio  Calculations*  To  predict  the  proton  affinities  and 
reaction  enthalpies  involved  in  this  study,  ab  initio  calculations 
were  carried  out  on  the  Gaussian  03  suile  of  programs46 
Molecular  geometries  were  preoptimized  ai  the  Hartree— Fock 
3-2  JG  level  of  theory,  and  vibraiional  frequencies  were 
calculated  to  confirm  thai  a  (local)  energy  minimum  had  been 
reached.  Subsequent  optimization,  energy,  and  vibrational 
frequency  calculations  were  performed  at  Hartree -Fock 
6-31+G(d,p)  level  or  density  functional  theory  (DFT)  level  with 
the  hybrid  functional,  I13LYP  6-31+G(d,p),  and  these  vibra¬ 
tional  frequencies  confirmed  that  the  optimized  structures  were 
minima  on  the  potential  energy  surface.  The  HF  zero-point 
vibrational  energies  (ZPVE)  were  scaled  by  0.9153  and  ihe  DFT 
zero-point  energies  were  scaled  by  0.9806. 47  ITie  corrected  zero- 
point  energies  were  added  to  ihe  stationary  point  energies.  No 
thermal  corrections  were  made  to  determine  the  reaction 
enthalpies  at  298  K,  and  the  values  reported  here  are  at  0  K. 
Reaction  enthalpies  at  0  K  were  calculated  by  subtracting  the 
sum  of  the  heats  of  formation  of  the  reactants  from  the  sum  of 
the  heats  of  formation  of  the  products.  Table  1  shows  the 
predicted  proton  affinities  and  reaction  enthalpies  at  0  K  of  the 
various  reaction  steps  considered  here.  A  more  detailed  ab  itiitio 
study  will  be  the  subject  of  a  future  publication. 

Results 

FUR  absorption  spectra  of  the  gas-phase  species  evolved 
when  ILs  react  with  white  fuming  nitric  acid  can  he  seen  in 
Figures  3-6.  In  Figure  3,  product  peaks  in  the  reaction  of 
l-butyl-3-methyl-imidazolium  azide  with  WFNA,  which  does 
not  ignite,  arc  shown.  To  improve  the  si  goal -to -noise  ratio,  an 
average  of  11  scans  (obtained  in  the  lime  range  270—620  ms) 
was  taken  after  the  spectrometer  was  triggered.  One  distinct 
NtO  and  two  HN3  peaks  are  delected  with  origins  at  2223.5, 
2139.6,  and  M49.4  cm"1,  respectively.  The  expected  HN3 
absorption  peak  ai  —3336  cm” 1  is  too  weak  to  be  discerned 
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Figure  3,  IK  product  spectrum  of  l-butyl-3-methyl-imidazoliunj 
{BMIM)  azide  reacting  v^dth  WFNA. 
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Figure  4.  IR  product  ftpectrutn  of  BMIM  dieyanamide  reacting  with 
WFNA. 
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Figure  5,  IR  product  spectrum  of  BMIM  dieyanamide  reacting  with 
WFNA.  The  evolution  of  CO;,  NzO<  and  UNCO  is  shown  here. 

from  ihe  recorded  spectral  noise.4*  All  other  peaks  in  the 
spectrum  have  been  assigned  to  gaseous  HNO3  and  NO;.  When 
l-butyl-3-methyl-irnidazofium  dieyanamide  is  reacted  with 
WFNA,  both  CO2  (origin  =  2348.4  cm”1)  and  N>0  (origin  — 
2222,6  cm  l)  are  evolved  before  ignition,  as  seen  in  Figure  4. 
Upon  closer  inspect  ion.  a  peak  is  seen  which  appears  in  between 
the  CO2  and  N2O  peaks  (Figure  5).  Figure  6  shows  the  result 
of  subtracting  the  NiO  and  C02  contributions  to  the  spectrum 
in  Figure  5.  after  having  normalized  the  N:G  peak  data  with 
Figure  3  and  the  C03  peak  with  an  ambient  spectrum.  The 
residual  absorption  feature  has  an  origin  at  approximately  2269 
cm  1  and  is  identified  as  isoeyanic  acid  (UNCO),  Absent  in 
the  spectra  of  the  preignition  vapor  phase  arc  any  other  gaseous 
combustion  products;  IICN.  1  lONO,  NO,  and  GO.  Absorptions 
by  lliO,  which  is  possible  according  to  our  mechanism  as 
discussed  later,  can  be  observed  in  Figure  4  near  —1600  and 
-3750  cm”1. 

The  time  profile  of  ihe  isoeyanic  acid  product  tracks  with 
the  NiQ  and  CO>  products,  as  seen  in  Figure  7,  Upon  systematic 
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substitution  of  eyano  groups  by  nitre  groups  on  the  d  icy  an  amide 
anion,  comp  orison  of  (be  reach  vily  of  sodium  dicy  an  amide 
(NaH(CN)2)T  silver  nitrocyanainide  (AgN(CN)NO^)T  and  potas¬ 
sium  dinitramide  (KN(N02)>)  with  WFNA  was  carried  out. 
Sodium  dicyanamide  reacts  violently  with  WFNA,  producing 
C XT,  N7O,  and  I  INCH  (see  Figure  ft)  before  igniting  (113 
^825  ms),  as  did  ail  the  heterocyclic  dicyanamide  ILs  with  dieir 
own  characteristic  ID  times  which  we  report  elsewhere.49  The 
observed  gas-phase  species  lime  profiles  in  the  reaction  of 
NaN(CNh  yah  and  WFNA  were  similar  to  that  in  Figure  7. 
The  nitrocyanamide  salt  does  not  ignite  and  FTIR  measurements 
showed  that  the  nitro  cyan  amide  anion  did  nor  react  to  produce 
any  detectable  gaseous  products  (in  fact,  only  the  initial  HNO;i 
and  NO:)  constituents  of  the  WFNA  arc  observed).  The 
dinitramidc  salt  reacted  slowly  without  ignition  to  evolve  N?G. 

In  comparing  the  spectra  of  all  dicyan  amides,  only  fire  CG?, 
N-rO,  and  iso  cyanic  acid  are  common  products.  Tire  observed 
CO^/N^O  product  signal  ratio  in  IL-dicyanamide  reactions  with 
WFNA  were  similar  to  that  in  sodium  dicyanamide's  reaction, 
suggesting  minimal  influence  of  the  cation  on  tire  reaction 
temperature  for  which  the  data  was  collected.  However,  for  both 
I -propargyl-.Vinelhyl-iTiLiday.ohum  dicynri  amide  arid  l-allyI-3- 
iirelhyTimidazolium  dicyanamide,  a  new  peak  emerges  which 
is  centered  at  1 3 00  cm-1  (Figure  9).  The  appearance  of  tills 
species  occurs  only  just  before  ignition  (within  15  to  45  msec) 
and  could  possibly  be  a  fragment  product  coming  from  reaction 
at  die  unsaturated  functional  group  of  die  cation. 

The  main  findings  of  the  biuret  reaction  test  were  as  follows. 
When  the  read  ion  of  NaNff.N)?  (90  mg)  with  34  wt  %  1  INO3 
(1  drop)  was  quenched  after  a  minute  rising  excess  2%  KOI  I 
solution,  a  cloudy  deep  purple-pink  complex  resulted  on  the 
addition  of  2—4  drops  of  \%  OuSO^  solution,  which  turned 
cloudy  pink  within  a  lew  minutes  and  remained  stable  for  many 
hours.  A  similar  result  was  also  obtained  for  1  -butyl- 3-mcthyl- 
imidazolium  dicyanamide/34  wt  %  HNO3  system.  Blank  runs 
in  which  no  acid  was  added  but  the  same  amount  of  KOI  ]  was 
used  produced  no  change  in  the  coloration  upon  Cu2+  addition. 
These  solutions  remained  light  blue.  The  pink  color  change  for 
the  dicyanamide/HNOj  systems  indicates  a  positive  biuret  test 
in  which  the  reaction  intermediate  (which  has  a  biuret  or  peptide- 
like  functionality)  reacts  with  the  Cu(II)  ions  to  form  a  Cu (In¬ 
coordination  complex.  When  the  same  tests  were  performed  for 
the  [  INO3  reaction  wiLh  l-methyl-d-amino-l  ,2,4,-lriav.olium 
nitrocyanamide  and  Tally  1-4-amino-l, 2,4,-  triazolium  nicrocy- 
anamide,  no  pink  coloration  was  observed. 

Discussion 

All  of  the  azides  when  treated  with  WFNA  or  inhibited  red 
fuming  nitric  acid  (TRFNA)  failed  to  ignite.  Displacement  by 
the  weaker  acid  on  excess  addition  of  nitric  acid  1  earls  to 
evolution  of  LlNn  into  the  gas  phase.  The  observed  N7O  (with 
coproduct  N?)  arises  as  a  result  of  decomposition  of  the  NON3 
intermediate^1  fonried  in  the  condensed- phase  interaction  of  N5” 
with  NOT  Hie  source  of  NO"1"  is  the  aiilokmi/atiou  reaction, 
(N?G,| NO+  4-  NO*-)  which  is  possible  since  small  amounts 
of  NO-2  can  be  present  in  WFNA  and  up  to  ^14%  is  dissolved 
in  WFNA  to  form  I  he  IRFNA.  7  he  initial  heal  release  from  I  he 
N3  ion  decomposition  to  form  N?  and  N?0  is  not  sufficient  in 
our  tests  to  further  decompose  the  heterocyclic  countercation 
to  canse  ignition,  even  when  (partial)  oxidation  of  the  cation 
may  concurrently  coTUribule  to  the  heal  release.44  CAUTION: 
The  product,  hy  dr  azoic  acid  (HN^)  is  a  highly  toxic  pungent 
smelling  compound  whose  vapor  can  cause  violent  headaches. 

!  INj  acts  as  a  non  cumulative  poison.  A  low  but  potentially  fatal 


Chambreau  el  al. 

dosage  of  NO?  by  inhalation  can  lead  to  lung  edema.  Low 
concentration  (4  ppm)  exposures  can  anesthetize  the  nose,  thus 
creating  a  potential  for  overexposure,  lr  must  also  be  emphasized 
thal  WFNA  and  IRFNA  are  violently  reactive,  hazardous,  and 
toxic  materials  and  only  small  quantities  should  be  used  under 
well- ventilated  conditions  with  cane  fully  planned  safety  proto¬ 
cols  and  protective  equipment  in  place.  The  present  experimental 
setup  of  Figure  1  (details  110L  shown)  allowed  for  die  proper 
containment  and  disposal  of  all  the  materials  Its,’ ted  and  their 
reaction  product. 

The  identification  of  IINCO  by  FTIR  and  the  elimination  of 
competing  possibilities  is  now  addressed.  Because  die  peak  at 
™22S9  cm“L  is  also  present  when  sodium  dicyanamide  salt  is 
reacted  with  WFNA,  it  cannot  be  a  product  resulting  from 
reaction  with  the  imidazolium  cation.  Tills  is  further  evidenced 
by  the  lack  of  high  frequency  hydrogen  stretching  modes 
expected  from  hydrocarbon  products  of  reaction  with  the 
imidazolium  cation,  and  which  are  also  ahsent  for  the  azide  IT. 
with  I  he  same  cation  (Figure  3),  Known  cyanogen 
(N=C— C— N)51  absorption  features  could  not  be  matched  with 
the  data  of  Figure  6,  and  the  presence  of  cyan  amide 
(I  frNO—N)52  is  also  ruled  out  based  on  the  fact  that  there  was 
no  evidence  for  absorption  al  1955  cm  1  due  Lo  the  tq(N  CN) 
symmetric  stretching  modeA^  Instead,  a  match  of  the  data  was 
found  with  the  spectrum  for  isoeyainc  acid  (HNCO). 54,55 
According  to  Ilerzberg  and  Reid*1  die  t^(N — C— O)  asymmetric 
Stic  telling  mode  centered  at  ■‘-’2274  cm-1  (and  more  recently 
reported  by  Steiner  and  co-workers  to  be  at  2268.9  cm  ~L  56  and 
Fisher  and  co- workers  to  be  at  2259  cm-157)  can  be  identified 
here  to  be  responsible  for  the  observed  strong  absorption.  HNCO 
also  has  two  smaller  peaks  centered  at  ^797  and  ^353 1  cm-1 
with  31%  and  24%  of  the.  peak  height  at  2274  cm~F  Even  with 
poor  signal -to- noise  ratio  in  the  800  cirT1  region,  the  former 
peak  for  the  fq(H— N— CG)  bend  is  just  discernible  in  Figure 
4,  while  the  latter  peak  for  ?q(H— N CO)  stretch  remains 
apparently  concealed  by  the  large  HNO3  peak  ar  around  3550 
cm  K  The  other  three  bands,  7/3,  and  respectively  at  1322, 
577,  and  656  cm  1  are  loo  weak,  and  in  any  case,  fire  former  is 
masked  hy  the  strong  I  ENO3  peaks  while  the  latter  S.wo  are  in 
I  he  detection  cutoff  vicinity  of  llie  present  detector.  Also,  the 
observed  2269  cm-1  feature  is  not  due  to  Lhe  ui  stretching 
features  of  fulminic  acid  (HCNO).  which  is  reported  to  be  at 
■~-j2196  cm-1  (and  thus  to  the  red  of  the  N20  absorption),5®^ 
or  of  cy  anic  acid  (HQCN)  reported  to  be  at  ’'-•2302  cm” 160  since 
its  corresponding  i^fHO— CN)  stretching  absorption  at  -'-1082 
cm-1  is  not  seen  in  our  data.61  Both  HCNO  and  HOCN  arc 
higher  energy  isomers  of  HNCO  and  are  not  expected  to  be  in 
any  significant  amounts  in  the  gas  phase62 

It  is  well-known  that  nitriles  can  undergo  acul  catalyzed 
hydration  reactions6^  because  the  protonaJ.ed  form  is  much  more 
receptive  Lo  nucleophiles  at  Lire  nitrile  carbon,  water  adds  slowly 
lo  the  carbon  and  then  rearranges  to  form  an  amide.  The  acid- 
calalyzed  reaction  with  Lire  dicyanamide  anion  can  occur  in  a 
similar  fashion,  but  when  reacting  with  WFNA,  fire  nucleophile 
is  NOT  instead  of  H?0,  and  die  nitrile  functional  group  is 
subsequently  converted  to  a  C(~ 0)NHNO?  group.  Because 
NO3-  is  a  better  nucleophile  than  water,64  the  rate  of  addition 
of  Nt>F  to  die  carbon  center  is  expected  to  be  faster  than  for 
water  in  WFNA  conditions. 

The  details  of  the  proposed  mechanism  for  the  reaction  of 
the  dicyanamide  anion  with  nitric  acid  are  presented  here.  The 
ItrsL  step  is  lo  protonate  the  dicyanamide  anion  at  the  nitrile 
nitrogen,  followed  hy  NO3  attack  al  the  electrophilic  carbon 
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and  NO;  migration  to  the  terminal  nitrogen  (sec  Tabic  1  for 
calculated  energetics): 


The  calculated  gas-phase  proton  affinity  of  the  nitrile  nitrogen 
is  about  7.6  kcal  mol  J  greater  than  that  for  the  central  nitrogen 
at  the  B3-LY P/5-3 1+G(d,p)  level  of  theory.  This  is  consistent 
with  previous  experimental  observations  for  the  preferred 
proto  nation  at  the  nitrile  nitrogen  in  the  topochemicaJ  solid- 
state  transformation  of  ammonium  dieyanamide  into  dicyan- 
di amide. 63  The  reduced  electron  density  at  the  imine-carbon 
makes  it  susceptible  to  nucleophilic  attack  by  NOC,  which  is 
calculated  to  be  exothermic  by  5.4  kcal  mol  l: 
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Subsequent  1 3- shift  of  NO;  to  the  terminal  nitrogen,  can  lead 
to  the  formation  of  a  C=0  double  bond  and  reaction  3  is 
calculated  to  be  exothermic  by  34.5  kcal  mol  "1: 
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energetically  favorable  as  are  many  of  the  proton  a  lion  steps. 
Other,  alternate  and  energetically  possible  mechanistic  routes 
for  DNB  formation*  particularly  in  nitric  acid  solutions,  could 
involve  either  attack  by  the  weaker  H;0  nucleophile  at  the 
imine-carbon  of  NCN=C=NH  followed  by  NO; 1  addition  and 
subsequent  deprotonation  to  form  nitrobiuret  (NB)*66  or  the 
NCN=C=NH  intermediate  may  first  undergo  another  protona- 
lion  followed  by  NOC  addilion/rearrangement  reactions  to  yield 
the  Nil.  These  possibilities  and  the  reaction  dynamics  of  DNB 
formation  and  decomposition  are  currently  being  studied  using 
ab  initio  methods  to  determine  reaction  enthalpies,  barrier 
heights,  transition  slate  properties,  etc. 

The  thermal  decomposition  of  DNB  has  previously  been 
shown  to  produce  CO;,  N;0,  and  I INCO  67  the  same  products 
that  we  have  observed  in  this  work: 
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The  above  reaction  sequence  is  similarly  possible  starting  at 
the  other  nitriJe  nitrogen  to  form  the  di nitrobiuret  (DNB)  anion: 
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Because  of  the  number  of  atoms  in  these  molecules,  the 
thermochemistry  was  not  calculated  tor  these  species.  However, 
the  enthalpies  of  reactions  4“6  are  likely  to  be  comparable  to 
reactions  1—3,  Final  proto  nation  at  the  central  nitrogen  will 
produce  DNB: 

The  results  of  the  present  calculations  tabulated  in  fable  1  show 
that  the  crucial  N03  addition  step  of  reaction  3  and  subsequent 
rearrangement  to  form  the  amide  linkage  of  reaction  4  are 


IOC  ^}| 


{1  (1 


H>"  NCO  +  HjN - Vlh  {£> 


UjO  -  N*0 


HNCO  +  Np  +•  CO: 


The  production  of  CO;,  N;G,  and  HNCO  from  the  derivative 
intermediate  and  N;0  and  ICO  from  the  nilramide  intermediate 
in  reaction  X  must  occur  through  a  multtstep  process,  and  the 
decomposition  of  these  chemically  labile  species  may  be  acid 
catalyzed.  The  proposed  dicyan amide/WFN A  reaction  mecha¬ 
nism  is  consistent  with  the  lack  of  any  1ICN  or  1BNCN  in  our 
gas-phase  FTIR  spectra.  Observation  of  only  isocyanic  acid, 
IINCX).  and  not  of  its  higher  energy  isomers  is  consistent  with 
their  known  thermal  instabilities  above  —100  °C68  Treatment 
of(NaNfCNh)  with  34%  and  68%  HN03  solutions  both  resulted 
in  similar  product  evolution  as  seen  in  Figure  5,  which  is 
consistent  with  DNB  formation  under  concentrated  acid  condi¬ 
tions,  Our  positive  biuret  reaction  tests  provide  strong  evidence 
for  a  peptide  like  functionality*  —  HN(0=C)NHC(=0)NH— f 
in  the  reaction  intermediate.  This  supports  the  stance  that  it  is 
the  reaction  of  dieyanamide  anion  that  initiates  the  decomposi¬ 
tion  of  the  tL  salt,  and  oxidation  of  the  cation  becomes  important 
only  later  after  ignition  has  occurred  in  WFNA. 

HNCO  is  known  to  react  in  the  gas  phase  with  water  to  form 
CO;  and  ammonia:'58'69 

H20  +  HNCO— NH3+C02  (9) 

Ammonia  can  react  very  rapidly  with  the  nitric  acid  to  form 
NH4NO3  which  could  be  why  it  was  not  detected  in  our  FTIR 
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Figure  fi.  Difference  product  spectrum  of  BM1M  dicy unamide  reacting 
with  WIN  A.  The  dark  line  is  the  spectrum  of  HNCO  when  the  C02 
and  N2()  contributions  are  subtracted  out. 


Figure  7.  Product  lime  profiles  in  the  react  ion  of  BMI M  dicyanamide 
with  WFNA:  CO;  (soJid).  HNCO  (dashed),  and  NjO  (dotted).  Data 
shown  is  for  case  when  mixture  did  not  ignite. 


cm' 

Figure  8,  IR  product  spectrum  of  sodium  dicyanamide  reacting  with 
WFNA,  Note  that  the  evolution  of  C02.  NT),  and  HNCO  is  similar  to 
Figure  5. 

setup.  Here,  corresponding  hydrolytic  decomposition  in  the 
condensed  phase,  H^O+  4*  UNCO  — ■  N1 14_  +  COj,  would  also 
be  consistent  with  lack  of  Nlh  gas-phase  signals.  Thermal 
decomposition  of  NII4NO3  is  known  to  be  accelerated  in  strong 
mineral  acids  and  can  lead  to  additional  NT)  formation.70  Thus, 
in  the  case  if  1 1NCO  reacts  to  completion  in  nitric  acid  solution, 
the  maximum  |C0?]/lN?0]  ratio  possible  would  be  2/3  as  shown 
below: 

NaN{CN)2+4HN03 — 2C02+3N2Q  +  2H:0  +  NaNO:> 

(10) 

Gas-phase  detection  of  l  INCO  in  our  drop-test  experiments 
would  suggest  lhat  the  observed  CO2/N1O  signal  represents  a 
ratio  less  than  2/3  during  the  preignition  phase  and  that  possibly 
further  reactions  of  HNCO  and/or  other  vapors  from  DNB 
decomposition  lead  to  lire  gas- phase  ignition  of  the  II,  fuels 
and  salts  studied  here  as  seen  in  Figure  2c.67’71  The  difference 
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Figure  9,  IR  product  spectrum  of  1-propargyl -3- methyl  imidazolium 
dicyari  amide  (solid)  and  UaUyF3-methyI-iitudazohuiti  die  van  amide 
(dashed)  reacting  with  WFNA  showing  possible  un saturated  species 
peak  near  1800  cm  K 


Figure  tO.  IR  product  spectrum  of  silver  cyanale  (AgOCN)  reacting 
with  WFNA.  Note  that  the  evolution  of  OO2  and  HNCO  is  similar  to 
lhat  shown  in  Figure  5,  hut  N20  production  is  less.  Ihe  dark  line  is 
the  difference  spectrum  for  HNCO, 

spectrum  in  Figure  6  shows  the  infrared  spectrum  of  HNCO 
formed  in  the  reaction  of  WIN  A  with  L -butyl- 3-methyl- 
imidazoltum  dicyanamide,  Figure  8  shows  product  spectrum  for 
WFNA  with  sodium  dicyanamide,  and  Figure  10  shows  the 
difference  spectrum  when  silver  cyanate  (AgOCN)  is  reacted 
with  WFNA.  The  ratio  of  intensities  of  CttyNT)  in  Figure  10 
is  different  from  that  in  Figure  8  indicating  a  different  reaction 
stoichiometry  than  for  the  dieyan amides,  with  less  N2G  forma¬ 
tion  since  its  source  can  only  be  from  hydrolysis  of  HNCO 
formed  in  the  exothermic  condensed  phase  reaction  of  AgOCN 
with  nitric  acid.  Observation  of  only  HNCO  in  the  gas  phase 
would  be  consistent  with  its  known  stability  relative  to  HGCN 
or  HCNO.  The  band  shape  at  2271  cm  1  in  Figure  ID  is  similar 
to  lhat  in  Figure  8.  further  confirming  that  indeed  UNCO  is 
released  in  the  dicyanamide/mtric  acid  systems.  No  visible 
llames  were  observed  in  AgOCN/WFNA  tests. 

'Hie  effect  of  nitro-substitution  for  eyano  groups  on  the 
dicyanamide  anion  reactivity  with  WFNA  was  investigated  by 
reacting  WIN  A  with  silver  nitrocyan  amide  or  potassium 
dinilramide.  Significant  gaseous  products  were  observed  for  the 
dinitramide  sail,  in  which  N*0  production  results  from  acid 
catalyzed  (and/or  thermal)  decomposition  of  the  dinitramic  acid72 
without  ignition: 

KN(N02)2+HN03  - 1  HN(N02)2  or  N(N(X)NOOH  |  + 

KNO^  (1 1} 

N{  NOONOOH  +  H + — N:0  +  H:0  +  NO:  (12) 

1 1,0  +  N02  —  ir  +  l  INQ3  ( 1 3) 
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HNCNCy  2  +  A — N?0  -h  HNO?  ( 14) 

The  ab  ini  hi)  calcul  Elions  indicate  l.hal  both  protonatiOTi  of  I  he 
centra!  nitrogen  followed  by  a  second  protonation  on  the  ni  tro¬ 
ts  xy  gen  or  double  protonation  on  the  nitre- oxygen  can  cause 
lhe  struct ure  to  dissociate  to  N?0,  IIjQ,  and  NO?  P  AH  for 
this  reaction  is  —205.5  keat  mol-1. 

HN(N02)2+H+  'NjO+HjO  I  NOj  (15) 

Silver  nitro cyanamide,  l-Tnethy1-4^amino-l,2,4-rria7oliuTn  ni- 
Lrocy an  ai  t  i  ide,  an  d  !  -al  ly  l  -4-ar  1 1  im )  - 1 , 2,4-  tri  a/.o  1  i  ui  1 1  t  i  i  In ) ey  ana- 
mide  failed  Lt>  ignite  in  our  acid  drop  tests,  From  Table  lt  Ihe 
lowest  energy  gas-phase  structure  for  the  protonation  of  the 
nitrocyanamide  anion  is  on  the  nitre-oxygen.  If  this  species  is 
subsequently  protonated  at  Lhe  central  nitrogen,  analogous  to 
die  dinitramidc  decomposition,  tlic  products  would  be  H2Ot  N^O, 
and  CN+: 

NCN(N02H)  +  H 1  —  N,0  +  H,0  +  CN (16) 

However,  the  energetics  of  forming  CN+  versus  NO?"1"  is  much 
less  favorable,  and  reaction  16  is  higher  in  energy  than  reaction 
15  by  at  least  210  kcal  mol-1,  which  is  likely  why  there  is  no 
N2()  evolution  upon  read  ion  of  TINO3  with  nilrocyanaTiLide.  In 
previous  work  on  decomposition  of  nilrocyanamide  in  strong 
mineral  acids,**  no  N2O  or  any  other  gaseous  products  were 
seen,  consistent  with  our  present  observations,  "Ibis  suggests 
that  unlike  in  the  above  din i Irani ide  case,  where  double 
pro  ton  at  ion  can  lead  to  exothermic  decomposition  to  give  N?0 
via  reaction  15,  here  this  step  docs  not  compete  with  protonation 
af  ihe  -UN  site.  Previously  Lhe  mechanism  for  ni  Irony  an  air  tide/ 
acid  reaction  in  which  protonalion  at  the  nitrile  nitrogen  occurs 
followed  by  nucleophilic  attack  by  H20  at  the  carbon  and 
subsequent  OH  migration  to  form  nirrourea  (NU)  has  been 
proposed.66  Nitro  urea  would  remain  in  Ihe  condensed  phase. 
Perhaps  the  amount  of  heat  generated  upon  formation  of 
nitrourea  is  insufficient  to  promote  higher  energy  reactions  and 
cause  either  gas  evolution  (i.e.T  N2G)  or  ignition.  Also,  Ihe 
decomposition  of  NU  or  further  nitration  of  NU  to  give 
dinitrourea  (DNU)  must  both  be  slow  in  our  drop  tests,  since 
GNU  is  thermally  unstable  and  should  also  give  gaseous  N?() 
and  CQ> 

The  allyk  prapargyl-,  and  (3-butcnyl)-  functionalized  mctliyl- 
imidazolium  dicyanamides  display  very  violent  ignitions  with 
WFNA.  giving  oil  a  white  Hash,  rather  than  a  typically  yellow 
Hash,  indicating  a  ho  Her  temperature  in  Lhe  ignition.  Evidence 
of  pyrolytic  decomposition  of  the  cation  was  also  seen  in  the 
high-speed  videos  in  which  copious  production  of  sootlike 
matter  upon  ignition  was  detected.  Initially,  these  functional 
groups  were  introduced  to  try  and  trigger  the  fuel-rich  cation 
to  ignite,  I  low  ever,  now  that,  the  clicyanm  11  ide  anion  is  lhe  trigger 
ibr  the  ignition,  the  cation  can  be  tailored  for  hie  best  overall 
properties  of  the  ionic  liquid.  For  example,  the  allyk  proparayk 
and  (3-huteny])-  functionalized  merhyl-imidazolmm  dicyana¬ 
mides  are  liquids  aJ  room  temperature  (rt),  while  the  (2- 
pentynyl)-  analogue  is  a  solid  at  rt.  The  anion  reacts  to  give  off 
significant  heat  in  the  fast  formation  of  C02  and  N^O,  which 
causes  Hubsequenl  decomposition  of  lhe  cation.  The  appearance 
of  an  unsaturated  C— C  stretch  might  indicate  that  the  addition 
of  unsaturated  functional  groups  causes  the  production  of 
smaller,  unsaturated  hydrocarbons  in  the  gas  phase  which  can 
greatly  aid  in  combustion  once  the  lemperalure  threshold  is 
reached  ("evidence  for  this  was  seen  as  secondary  ignition  flashes 
in  the  high-speed  video  recordings).  Tins  peak  is  only  seen  in 
the  ally  I-  and  propargyl- functionalized  melhyl-imidaz.olium 
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d icy anamide  species,  and  not  m  the  corresponding  (3-butcnyl)- 
system,  and  thus  more  work  is  needed  l.o  uncover  the  details  of 
how  substitution  affects  die  decomposi  Liou/oxidatiou  mecha¬ 
nisms  of  the  heterocyclic  cation. 

Conclusions 

The  hypcrgoLic  ignition  of  dicyanamide  ionic  liquids  upon 
reaction  with  fuming  nitric  acid  has  been  demonstrated,  and  it 
has  been  shown  that  the  anion  reaction  is  responsible  for 
initiating  the  ignition.  The  evolution  of  C02,  NsO,  and  HNCO 
during  prei  guidon  indicates  a  complex  reaction  median  ism, 
perhaps  Lh  rough  dinitrobiuret  and  nil  rani  ide  inlermediales. 
Evidence  for  die  former  was  con  dimed  in  the  biuret  reaction 
test  which  was  only  positive  for  the  dicyanamide  reactions  and 
not  with  any  of  Lhe  other  anions  studied  here.  In  Lhis  work  we 
have  shown  that  the  initial  reactivity  of  the  anions,  -N(CN)21 
“NNO-zfCN),  and  ”N(NOo)2r  with  WFNA  varies  dramatically 
with  Nt>2  substitution.  ■Similarly,  it  may  also  be  possible  to  affect 
Lhe  cation  reaclivily  by  suitable  chemical  substitutions,  so  as  to 
tailor  its  oxidation  rate.  Further  spectroscopic  probing  of 
products  and  intermediates  during  the  preignition  and  ignition 
phases  mid  additional  ab  iniLio  calculations  will  be  valuable  in 
providing  mechanistic  insight  on  how  IL- dicyanamides  undergo 
combustion  when  treated  with  WFNA  and  IRFNA. 
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Combined  data  of  photoelectron  spectra  and  photo  ionization  efficiency  curves  in  the  near  threshold  ionization 
region  of  isolated  ion  pairs  from  [emiirOtTtfcN],  [eimmJfPfiN],  anti  [diMpim][TfaN]  ionic  liquid  vapors  reveai 
small  shifts  in  Ihe  ionization  energies  of  ion-pair  systems  due  to  cation  ;*nd  anion  substitutions.  Shifts  toward 
higher  binding  energy  following  anion  substitution  are  attributed  to  increased  electronegativity  of  the  anion 
itself,  whereas  sliifts  toward  low  er  binding  energies  following  cation  substitution  are  attributed  to  ail  increase 
in  the  cation— anion  distance  that  causes  a  !ow?er  Coujombic  binding  potential.  The  predominant  ionization 
mechanism  in  (he  near  threshold  photon  energy  region  is  identified  as  dissociative  ionization,  involving  the 
dissociation  of  the  ion  pair  and  the  production  of  intact  cations  as  the  positively  charged  products. 


Introduction 

RoonHeisiperatUTEi  ionic  liquids  (Ehs)  are  sall.s  with  an  telling 
point  near  or  below  ambient  temperatures,1  They  are  typically 
based  on  bulky  mulliatom  cations  anti  anions  dial  resirici  the 
cation— anion  distance  of  closest  approach,  thus  limiting  their 
Con  torn  hie  attraction  arid  consequently  also  lowering  their 
melting  temperature  even  below  room  temperature.  Melting 
temperatures  of  room  temperature  ILs  can  therefore  be  much 
lower  than  those  of  sails  based  oti  atomic  ions,  such  as  NaCl, 
which  has  a  melting  temperature  of  1074  K..  In  addition  to 
having  low  molting  temperatures,  ILs  have  a  number  of  unique 
properties,  including  excellent  solvation  ability,  extraordinarily 
low  volatility,  and  high  thermal  stability. 1-5  Eiecause  rhe  physical 
and  chemical,  properties  of  ILs  can  be  varied  by  the  choices  of 
the  cation  and  anion  pair,1  [hey  can  be  used  for  many 
applications  such  as  fuel  cells,3  s  batteries, 6,7  and  solar  cells.8-5 
Furthermore,  ILs  are  also  used  in  ionic  propulsion10,11  and  for 
hypergohe  fuels1  ?  and  various  chemical  extraction  applications.11 
Hie  vast  number  of  possible  anion— cation  combinations1''1  makes 
it  valuable  to  predict  IL  properties  rationally  on  the  basis  of 
rhei r  constituents,15  and  considerable  experimental  and  theoreti¬ 
cal  work  has  been  performed  to  predict  and  measure  IL 
properties. Ab  initio  calculations,  used  to  predicl  Ihe  sl.ru  dure 
and  interactions  of  ILs,  are  typically  performed  on  isolated  "gas- 
phase'"  systems.  ‘ 71 It  is  therefore  attractive  to  study  isolated 
ion  pairs  and  lo  examine  basic  J I ,  decomposition  pathways,  such 
as  dissociation  or  ionization.  Because  ILs  arc  generally  perceived 
as  “green  solvents",  their  decomposition  pathways  can  lead 
to  a  variety  of  products  lhal  need  lo  be  considered  in  specific 
applications.  Given  recenl  works-"  demonstrating  the  distillation 
of  selected  families  of  ILs,  we  subsequently  prepared  molecular 
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beams  of  IL  vapor,  providing  a  gas-phase  target  of  isolated  ion- 
pair  species  to  obtain  the  pboi.oelectron  (PE)  spec  trim  i  of 
[enihuJlTl^N]  vapor. J,j  The  direct  comparison  of  our  [eTiiiinJ- 
[TfzN]  vapor  PE  spectrum  with  previously  measured  PE  spectra 
of  a  [bmimlfTfigN]  liquid  surface7 1>7;'  provided  important  evi¬ 
dence  of  the  ion-pah  nature  of  IL  vapor,"0  Other  recent  studios 
aimed  al.  in vesti  gating  I  be  proper  lies  of  IL  vapor  also  include 
line  of  sight  mass  spectroscopy,2 2-24  IR  spectroscopy,23  and 
matrix- isolated  FTIR  spectroscopy*  1 6  Moreover,  mass  spectra 
of  cation-  and  anion-substituted  ILs20  showed  that  the  ionization 
of  these  isolated  IL  ion  pairs  by  absorption  of  a  single  23.2  eV 
photon,  which  is  >14  eV  above  the  ionization  threshold,  results 
in  fragmentation  with  a  predominant  product  of  the  intact 
callous.  Interestingly,  no  signal  is  delec  ted  for  the  ionized  parenl 
ion-pair  mass;  that  is,  Ihe  removal  of  an  electron  frotu  the  ion- 
pair  system  results  in  the  dissociation  of  the  ion  pair.  Because 
those  results  were  obtained  with  a  photon  energy  that  is  14.3 
cV  higher  than  the  energy  needed  to  remove  an  electron  from 
the  ion  pair,  which  could  provide  ample  energy  to  detach  an 
electron  as  wrell  dissociate  the  molecule,  near  threshold  inves¬ 
tigations  are  required.  I 'id d- induced  IL  ionization  studies 
recently  provided  evidence  of  the  possible  formation  of  undis¬ 
sociated  parent  ion  systems  with  the  combined  mass  of  the 
cation  and  anion.27 

Here  we  present  a  study  of  the  near  threshold  soft  photoion¬ 
ization  of  isolated  ion  pairs  with  tunable  photon  energies 
b  el  ween  8.3  and  9,5  eV.  Electronic  si  rue  Eure  changes  following 
cation  and  anion  substitutions  are  presented  and  discussed  in 
terms  of  the  observed  binding  energy  shifts.  Starting  from  ihe 
c  o  m  m  on  1  y  used  [  I  -  eth  y  1  -  3  -m  erh  y  1  im  idazo  I  ium  ]  [b  is  ( rri  fi  no  - 
ronielhyIsuliimyl)iinide)J  ([eiriiTnjjTLNj)  rootii- temperature  IL, 
shown  in  Figure  1,  the  [emiroj 1  cation  is  substituted  by  die 
bulkier  imidazolium- based  1,2 -dimethyl- 3- propylimidazolium 
ion  (Jdmpim]+),  and  die  [TEN]-  anion  is  substituted  by  die 
bisfpentafhiorocthylsulfonyLimide  ion  (TPf2Nl  h  which  lias  a 
higher  fluorine  content.  Single-photon  near  threshold  photoion¬ 
ization  mass-spectra  and  photoionization  efficacy  (PIE)  curves, 
obtained  at  the  Advanced  1  right  Source  chemical  dynamics  beam 
line  in  the  near  threshold  region,  are  presented  for  the  c  at  ion  - 
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Combined  data  of  photoelectron  spectra  and  photoionization  efficiency  curves  in  the  near  threshold  ionization 
region  of  isolated  ion  pairs  from  [emim][Tf2N],  [emim][Pf2N],  and  [dmpim][Tf2N]  ionic  liquid  vapors  reveal 
small  shifts  in  the  ionization  energies  of  ion-pair  systems  due  to  cation  and  anion  substitutions.  Shifts  toward 
higher  binding  energy  following  anion  substitution  are  attributed  to  increased  electronegativity  of  the  anion 
itself,  whereas  shifts  toward  lower  binding  energies  following  cation  substitution  are  attributed  to  an  increase 
in  the  cation— anion  distance  that  causes  a  lower  Coulombic  binding  potential.  The  predominant  ionization 
mechanism  in  the  near  threshold  photon  energy  region  is  identified  as  dissociative  ionization,  involving  the 
dissociation  of  the  ion  pair  and  the  production  of  intact  cations  as  the  positively  charged  products. 


Introduction 

Room-temperature  ionic  liquids  (ILs)  are  salts  with  a  melting 
point  near  or  below  ambient  temperatures.1  They  are  typically 
based  on  bulky  multiatom  cations  and  anions  that  restrict  the 
cation— anion  distance  of  closest  approach,  thus  limiting  their 
Coulombic  attraction  and  consequently  also  lowering  their 
melting  temperature  even  below  room  temperature.  Melting 
temperatures  of  room  temperature  ILs  can  therefore  be  much 
lower  than  those  of  salts  based  on  atomic  ions,  such  as  NaCl, 
which  has  a  melting  temperature  of  1074  K.  In  addition  to 
having  low  melting  temperatures,  ILs  have  a  number  of  unique 
properties,  including  excellent  solvation  ability,  extraordinarily 
low  volatility,  and  high  thermal  stability.1-9  Because  the  physical 
and  chemical  properties  of  ILs  can  be  varied  by  the  choices  of 
the  cation  and  anion  pair,1  they  can  be  used  for  many 
applications  such  as  fuel  cells,3-5  batteries,6,7  and  solar  cells.8,9 
Furthermore,  ILs  are  also  used  in  ionic  propulsion10,11  and  for 
hypergolic  fuels12  and  various  chemical  extraction  applications.13 
The  vast  number  of  possible  anion— cation  combinations14  makes 
it  valuable  to  predict  IL  properties  rationally  on  the  basis  of 
their  constituents,15  and  considerable  experimental  and  theoreti¬ 
cal  work  has  been  performed  to  predict  and  measure  IL 
properties.16  Ab  initio  calculations,  used  to  predict  the  structure 
and  interactions  of  ILs,  are  typically  performed  on  isolated  “gas- 
phase”  systems.17,18  It  is  therefore  attractive  to  study  isolated 
ion  pairs  and  to  examine  basic  IL  decomposition  pathways,  such 
as  dissociation  or  ionization.  Because  ILs  are  generally  perceived 
as  “green  solvents”,  their  decomposition  pathways  can  lead 
to  a  variety  of  products  that  need  to  be  considered  in  specific 
applications.  Given  recent  works19  demonstrating  the  distillation 
of  selected  families  of  ILs,  we  subsequently  prepared  molecular 
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beams  of  IL  vapor,  providing  a  gas-phase  target  of  isolated  ion- 
pair  species  to  obtain  the  photoelectron  (PE)  spectrum  of 
[emim][Tf2N]  vapor.20  The  direct  comparison  of  our  [emim]- 
[Tf2N]  vapor  PE  spectrum  with  previously  measured  PE  spectra 
of  a  [bmim][Tf2N]  liquid  surface21,22  provided  important  evi¬ 
dence  of  the  ion-pair  nature  of  IL  vapor.20  Other  recent  studies 
aimed  at  investigating  the  properties  of  IL  vapor  also  include 
line  of  sight  mass  spectroscopy,23,24  IR  spectroscopy,25  and 
matrix-isolated  FTIR  spectroscopy.26  Moreover,  mass  spectra 
of  cation-  and  anion-substituted  ILs20  showed  that  the  ionization 
of  these  isolated  IL  ion  pairs  by  absorption  of  a  single  23.2  eV 
photon,  which  is  >14  eV  above  the  ionization  threshold,  results 
in  fragmentation  with  a  predominant  product  of  the  intact 
cations.  Interestingly,  no  signal  is  detected  for  the  ionized  parent 
ion-pair  mass;  that  is,  the  removal  of  an  electron  from  the  ion- 
pair  system  results  in  the  dissociation  of  the  ion  pair.  Because 
those  results  were  obtained  with  a  photon  energy  that  is  14.3 
eV  higher  than  the  energy  needed  to  remove  an  electron  from 
the  ion  pair,  which  could  provide  ample  energy  to  detach  an 
electron  as  well  dissociate  the  molecule,  near  threshold  inves¬ 
tigations  are  required.  Field-induced  IL  ionization  studies 
recently  provided  evidence  of  the  possible  formation  of  undis¬ 
sociated  parent  ion  systems  with  the  combined  mass  of  the 
cation  and  anion.27 

Here  we  present  a  study  of  the  near  threshold  soft  photoion¬ 
ization  of  isolated  ion  pairs  with  tunable  photon  energies 
between  8.3  and  9.5  eV.  Electronic  structure  changes  following 
cation  and  anion  substitutions  are  presented  and  discussed  in 
terms  of  the  observed  binding  energy  shifts.  Starting  from  the 
commonly  used  [l-ethyl-3-methylimidazolium][bis(trifluo- 
romethylsulfonyl)imide)]  ([emim][Tf2N])  room-temperature  IL, 
shown  in  Figure  1,  the  [emim]+  cation  is  substituted  by  the 
bulkier  imidazolium-based  l,2-dimethyl-3-propylimidazolium 
ion  ([dmpim]+),  and  the  [Tf2N]-  anion  is  substituted  by  the 
bis(pentafluoroethylsulfonyl)imide  ion  ([Pf2N]-),  which  has  a 
higher  fluorine  content.  Single-photon  near  threshold  photoion¬ 
ization  mass-spectra  and  photoionization  efficiency  (PIE)  curves, 
obtained  at  the  Advanced  Light  Source  chemical  dynamics  beam 
line  in  the  near  threshold  region,  are  presented  for  the  cation- 
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Figure  3,  Photoeleetron  binding  energy  spectra  of  IL  vapor  taken  with 
23.2  eV  photons  (53.7  nm)  for  [emim][Tf?N],  [emim][Pf2N],  and 
[dtnpim][Tf2N]. 


Figure  4.  Detail  of  the  PE  spectra  near  the  ionization  threshold  of 
|emint][Tf;N],  femim ||Pf*N |,  and  [dmpim][Tf:N]  obtained  with  23.2 
eV  photons  (53.7  nm). 


[emim]  [TEN]  and  [emimHPt^N]  on  the  normalized  soak.  The 
solid  lines  in  Figure  4  show  fils  of  the  PE  spectra  in  the  neat 
threshold  region.  We  use  Ihe  simplest  threshold  function  to  lit 
the  data,  which  vanishes  below  the  threshold  energy  and  rises 
linearly,  as  cl(E  —  JSy,  for  binding  energies  beyond  the  threshold 
ionization  energy.  Here  £b  is  the  fitted  threshold  energy  and  a 
is  a  tilled  overall  normalization  factor,  which  depends  on  the 
photoionization  cross  section,  photon  flux,  and  sample  density. 
This  simple  threshold  function  is  convoluted  with  a  250  meV 
Gaussian.  The  width  of  the  Gaussian  was  also  fitted  to  the  data 
and  reflects  the  combined  effect  of  internal  excitation  of  the 
hot  ion  pairs,  the  bandwidth  of  the  light  source,  and  the 
experimental  response  function  of  the  PE  spectrometer.  We  find 
that  a  similar  trend  of  binding  enrgy  shifts  observed  for  peak 
positions  can  be  identified  in  the  fitted  threshold  ionization 
energies.  Whereas  for  femim][Tf2M],  the  threshold  ionization 
energy  is  found  to  be  8,9  ±  0.2  eV,  the  threshold  for 
lemimliPENJ  is  shifted  to  higher  binding  energies  9,0  ±  0,2 
eV,  and  for  [dmpim)[Tf2N|,  the  threshold  ionization  energy  is 
shifted  to  lower  binding  energies  of  8,7  ±  0.35  eV,  The  changes 
in  the  electronic  structure  due  to  the  cation  and  anion  substitu¬ 
tions  are  small  and  can  be  mainly  characterized  by  these  slight 
binding  energy  shifts.  The  main  sources  of  error  in  this 
measurement  are  the  limited  resolution  of  the  PE  spectrometer 
for  fast  electrons  emitted  near  the  ionization  threshold  as  well 
as  the  noise  due  to  background  electrons  in  this  PE  region.  The 
trends  of  the  small  binding  energy  shifts  suggested  by  the 
features  of  the  measured  PE  spectra  are  reaffirmed  in  l lie 


mass  (amu) 


Figure  5,  Near-threshold  photoionization  mass  spectra  collected  with 
photon  energies  in  the  K.6  to  9.5  eV  range  for  (a)  [emim][Tf:R|.  (b) 
[cmimltPfiNl.  and  (c)  [dmpimllTfzNl, 


following  by  the  higher  resolution  tunable- wavelength  PIE 
spectra  measurements. 

The  photoionizalion  of  the  IL  ion  pairs  with  23,2  eV  photons, 
which  is  used  to  obtain  the  PE  spectra  and  results  in  dissociative 
ionization  and  a  predominant  mass  peak  of  the  respective  intact 
cation  signal  in  photoionizalion  mass  spectra.20  With  the 
implementation  of  tunable  wavelengths  at  the  ALS  chemical 
dynamics  beamline,  it  is  possible  to  tune  the  ionizing  photon 
energies  near  the  thresholds  of  ionization  that  are  revealed  by 
the  PE  spectra.  Therefore,  a  relatively  soft  ionization  mechanism 
is  used  for  examining  the  stability  of  isolated  IL  ion  pairs.  Figure 
5  shows  near- threshold  photoionizalion  mass  spectra  of 
Lemimj[Tf\N|.  [emim] [PEN],  and  [dmpim][Tf2NJ  IL  vapors 
obtained  in  the  near-ionization  threshold  region  with  photon 
energies  between  8,6  and  9.5  eV,  Similar  to  mass  spectra 
obtained  with  higher  photon  energies,  we  find  that  the  predomi¬ 
nant  positively  charged  products  of  the  considered  ion- pair 
systems  are  the  respective  intact  cations,  that  is,  [emim]'  at 
mass  1 1 1  amu  and  [dmpim] '  at  mass  139  amu.  It  is  also  possible 
to  observe  the  corresponding  ]3C  isotope  peaks  at  112  and  140 
amu.  The  relative  yields  of  the  13C  isotope  peaks  are  in 
agreement  with  the  natural  isotope  abundance  for  |emim]+  with 
eight  carbon  atoms  and  [dmpim]+  with  ten  carbon  atoms. 
Despite  the  fact  that  soft  photo  ionization,  with  near- threshold 
photon  energies,  was  used,  no  parent  mass  peaks  were  observed 
(not  shown).  The  excellent  level  of  signal -to- noise  allows  us  to 
assess  an  upper  limit  of  belter  than  M:1G0G  for  the  relative 
yield  of  parent  ions  at  the  mass  of  the  undissociated  ion  pairs 
compared  with  the  observed  yield  of  dissociated  cation  mass. 
This  indicates  that  the  threshold  ionization  energies  of  these 
isolated  ion  pairs  coincide  with  the  appearance  energies  of  the 
intact  cations  produced  by  dissociative  ionization.  It  is  therefore 
concluded  that  the  removal  of  the  highest  occupied  molecular 
orbital  (HOMO)  electron  from  the  IL  ion  pairs  considered  here 
leads  to  dissociation  and  formation  of  the  intact  cation. 

To  provide  a  more  precise  determination  of  the  small  binding 
energy  shifts  at  the  threshold  ionization  energies,  for  the  cation 
and  anion  substitutions,  we  record  PIE  curves  as  a  function  of 
photon  energy  in  the  near- threshold  ionization  region.  On  the 
basis  of  the  fact  that  only  dissociative  ionization  is  observed  in 
the  ionization  threshold  region,  as  identified  by  PE  spectroscopy, 
we  assume  that  the  appearance  energy  of  the  intact  cation 
coincides  with  the  ionization  energy.  Therefore,  PIE  curves  are 
constructed  by  recording  the  yield  of  the  respective  intact  cation 
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Figure  6.  Near- threshold  PIE  spectra,  showing  binding  energy  shifts 
due  to  cation  and  anion  substitutions. 

signal  normalized  to  the  ALS  photon  flux,  as  a  function  of 
photon  energy.  Figure  6  shows  the  near- threshold  PIE  curves, 
recorded  for  |emhnl[TEN],  [ending PEN],  and  Jdmpim][TEN], 
PIE  curves  can  be  compiled  with  the  integral  of  the  PE  spectra. 
Therefore,  we  lit  the  near  threshold  PIE  curves  with  a  function 
that  is  the  integral  of  the  linear  threshold  function  we  success¬ 
fully  used  for  fitting  the  near- threshold  PE  spectra,  that  is, 
vanishing  yield  below  the  threshold  energy  and  a  quadratic  rise 
of  the  ionization  efficiency  above  the  threshold  energy.  There¬ 
fore  the  functional  form  is  l/2a(E  —  £^>)2,  where  a  is  a  lilted 
overall  normalization  factor,  E  is  the  photon  energy,  and  ED  is 
the  filled  threshold  parameter.  Interestingly,  this  functional  form 
that  successfully  describes  our  data  is  also  known  as  the  Fowler 
law,  which  was  originally  developed  lo  descrihe  photoemission 
from  inetals?t  but  was  found  to  be  useful  for  other  systems  as 
well.32  Similarly  to  the  PE  spectrum  lilting  procedure,  this 
function  is  convoluted  with  a  Gaussian  resolution  function,  fitted 
to  be  narrower,  50  meV.  and  reflects  She  better  resolution  of 
the  monochromator  at  the  9. 0.2.3  beamline.  The  three  data  sets 
shown  in  Figure  6  are  shown  lo  be  normalized  by  their  lilted 
overall  normalization  factor,  a.  The  solid  lines  in  Figure  6  are 
die  fitted  threshold  functions,  corresponding  to  threshold  ioniza¬ 
tion  energies  of  8.72  ±  0.03,  8.78  ±  0.04,  and  8.59  +  0.03  eV 
for  [emim][TEN],  [emim][PEN],  and  [dmpim][TEN],  respec¬ 
tively.  The  fitted  threshold  ionization  energies  are  in  good 
agreement  with  the  analysis  of  the  PE  spectra  within  the  error 
bars  of  the  PE  spectra  analysis.  Although  within  the  experi¬ 
mental  PE  error  bars  of  0.2  eV,  the  ionization  energies  obtained 
by  the  PE  spectra  analysis  are  systematically  higher  than  PIE 
results,  which  is  possibly  due  to  a  small  shift  on  the  order  of 
0. 1  eV  in  the  PE  spectrometer  calibration.  The  results  reaffirm 
with  higher  precision  the  trends  of  the  small  binding  energy 
shifts  observed  in  the  PE  spectra.  It  is  important  to  note  that 
for  the  three  IL  experiments,  the  reservoir  of  the  ion  source  is 
kept  at  similar  temperatures  around  500  Kf  and  thus  the  observed 
shifts  are  not  expected  to  be  due  lo  different  molecular  beam 
temperatures.  The  substitution  of  the  [TEN]  anion  by  [PEN]  ' 
results  in  a  small  +60  meV  shift  toward  higher  electron  binding 
energies.  We  suggest  that  this  trend  may  be  understood  in  terms 
of  increased  electronegativity  due  to  the  substitution  of  the 
irifluoro methyl  by  pentafluoroethyl  groups,  leading  to  higher 
binding  energy  of  the  electron  to  the  anion.  In  the  case  of 
[emim] cation  substitution  by  [dmpimj  \  the  —130  meV  shift 
is  toward  lower  binding  energies.  The  addition  of  the  2-methyl 
group  replacing  a  hydrogen  atom  and  the  replacement  of  3-ethyl 
by  the  bulkier  3-propyl  group  can  act  to  alter  the  spacing 
between  the  ions,  increasing  the  cation— anion  distance  and 
therefore  reducing  the  Coulombic  interaction  between  the  cation 
and  the  electron  of  the  anion.  Therefore,  we  suggest  that  this 


Strasser  et  ill, 

particular  cation  substitution  can  lead  to  the  observed  binding 
energy  shill  toward  lower  binding  energies  of  the  anions 
electron. 

Conclusions 

Isolated  ion  pairs  in  the  IL  vapor  of  three  IL  species  were 
studied  with  a  tunable  synchrotron  light  source  in  the  near 
threshold  ionization  region,  between  8.3  and  9.5  eV  photon 
energies  as  well  as  with  23.2  eV  photons  from  a  “table-top” 
high- order  harmonic  generation  light  source.  The  combined  PE 
spectra,  near- threshold  mass  spectra,  and  PIE  curves  allow  ns 
to  determine  ionization  threshold  energies  and  to  detect  small 
binding  energy  shifts  due  to  cation  and  anion  substitutions. 
Theoretical  calculations  are  needed  lo  provide  further  insight 
into  the  electronic  structure  and  the  observed  binding  energy 
shifts  of  gas- phase  IL  ion  pairs. 

Photoionization  of  the  examined  IL  vapor  predominantly  leads 
to  dissociative  ionization.  Photodetachment  of  the  HOMO 
electron  results  in  the  dissociation  of  the  ionized  ion-pair  system 
into  a  neutral  fragment  and  the  intact  cation  that  is  detected  in 
the  mass  spectra  measurements.  Previous  works  using  pulsed 
field  ionization  of  similar  IL  species  have  showed  that  the 
ionized  [cation] [anion]  mass  peak  can  be  delected  in  positive 
ion  mass  spectra,27  However,  in  the  present  work,  even  with 
the  introduction  of  soft  photoionization,  there  is  no  evidence 
of  the  undissoejated  parent  system  in  the  recorded  mass  spectra. 
One  possibility  is  that  the  parent  system  formed  by  photoion- 
ization  is  unstable  and  spontaneously  dissociates  on  a  time  scale 
that  is  much  faster  than  the  typical  tens  of  microseconds  between 
the  quasi-continuous  ALS  light  source  and  the  pulsed  TOE  ion 
extraction  m  our  experiments. 
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The  heats  of  vaporization  of  the  room  temperature  ionic  liquids  (RTILs)  A-butyl-A-methylpyrrolidinium 
bistrifluorosulfonylimide,  A-butyl-A-methylpyrrolidinium  dicyanamide,  and  l-butyl-3-methylimidazolium 
dicyanamide  are  determined  using  a  heated  effusive  vapor  source  in  conjunction  with  single  photon  ionization 
by  a  tunable  vacuum  ultraviolet  synchrotron  source.  The  relative  gas  phase  ionic  liquid  vapor  densities  in  the 
effusive  beam  are  monitored  by  clearly  distinguished  dissociative  photoionization  processes  via  a  time-of- 
flight  mass  spectrometer  at  a  tunable  vacuum  ultraviolet  beamline  9. 0.2. 3  (Chemical  Dynamics  Beamline)  at 
the  Advanced  Light  Source  synchrotron  facility.  Resulting  in  relatively  few  assumptions,  through  the  analysis 
of  both  parent  cations  and  fragment  cations,  the  heat  of  vaporization  of  A-butyl-A-methylpyrrolidinium 
bistrifluorosulfonylimide  is  determined  to  be  A//vap(298.15  K)  =  195  db  19  kJ  mol-1.  The  observed  heats  of 
vaporization  of  l-butyl-3-methylimidazolium  dicyanamide  (A//vap(298.15  K)  =  174  db  12  kJ  mol-1)  and 
A-butyl-A-methylpyrrolidinium  dicyanamide  (A77vap(298.15  K)  =  171  db  12  kJ  mol-1)  are  consistent  with 
reported  experimental  values  using  electron  impact  ionization.  The  tunable  vacuum  ultraviolet  source  has 
enabled  accurate  measurement  of  photoion  appearance  energies.  These  appearance  energies  are  in  good 
agreement  with  MP2  calculations  for  dissociative  photoionization  of  the  ion  pair.  These  experimental  heats 
of  vaporization,  photoion  appearance  energies,  and  ab  initio  calculations  corroborate  vaporization  of  these 
RTILs  as  intact  cation— anion  pairs. 


Introduction 

Ionic  liquids  are  classified  as  ionic  salts  that  have  melting 
points  at  and  below  T  =  100  °C.  Thus,  room  temperature  ionic 
liquids  (RTILs)  are  liquids  under  ambient  conditions.  RTILs 
represent  a  new  class  of  solvents  that  potentially  can  replace 
highly  volatile  organic  solvents  and  hydrazine-based  hypergolic 
fuels.1-4  RTILs  typically  have  large,  multiatom  cations  such  as 
ammonium,  pyrrolidinium,  imidazolium,  triazolium,  and  tetra- 
zolium  ions  that  have  diffuse  positive  charge  distributions  and 
nonpolar  alkyl  functional  groups  (Figure  1).  The  diffuse  nature 
of  the  ionic  charges,  combined  with  the  large  effective 
cation— anion  distance  in  RTILs  results  in  relatively  low  lattice 
energies.  By  contrast,  a  typical  ionic  salt,  NaCl  (melting  point 
T=  801  °C),  has  very  localized  charge  distributions  on  each 
ion,  and  the  separation  between  cation  and  anion  is  small; 
therefore,  a  strong  Coulombic  attraction  exists  between  sodium 
and  chloride  ions  in  the  condensed  phase.  As  a  result,  the  lattice 
energy  of  NaCl  is  much  higher  than  that  for  RTILs,  and  the 
phase  transition  from  ionic  solid  to  liquid  occurs  at  a  much 
higher  temperature  in  NaCl  than  for  RTILs. 
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Figure  1.  A-butyl-A-methylpyrrolidinium  bistrifluorosulfonylimide 
[pyrri4]+[NTf2]_,  a  room  temperature  ionic  liquid. 


Due  to  their  ionic  nature,  RTILs  were  originally  thought  to 
have  essentially  no  vapor  pressure.  Recent  studies  of  RTILs  in 
a  vacuum  have  demonstrated  that  some  RTILs  can  be  distilled 
in  a  vacuum  with  little  thermal  degradation,5-8  as  noted  in  more 
detail  below.  Transpiration  and  Knudsen  methods9  have  also 
been  applied  to  vaporize  and  recover  RTIL  condensates.10-13 
To  evaluate  the  thermal  stability  of  a  compound,  studies  such 
as  differential  scanning  calorimetry  (DSC)  and  thermal  gravi¬ 
metric  analysis  (TGA)  are  typically  carried  out  at  atmospheric 
pressure.11,12,14-16  Recent  TGA  studies  have  shown  that  signifi¬ 
cant  mass  loss  of  RTILs  can  occur  at  temperatures  well  below 
the  onset  decomposition  temperatures  reported  in  DSC  studies.16 
Such  loss  in  mass  could  represent  depletion  through  vaporization 
of  intact  RTILs. 

The  ionic  nature  of  RTILs  may  lead  to  dramatically  different 
vaporization  mechanisms  than  for  molecular  liquids.  For  protic 
RTILs  (having  an  H  bound  to  a  nitrogen  in  the  cation  ring), 
typically  a  proton  is  transferred  to  the  anion  during  vaporization 
and  neutral  species  are  formed:7 
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[AH]+[B]-  —  A  +  BH  (1) 

These  neutrals  can  be  volatilized,  and  upon  condensation,  the 
RTIL  can  be  reformed  by  proton  exchange  to  the  weaker  acid 
(the  reverse  of  reaction  1  above).  The  groundbreaking  study 
by  Earle  et  al.5  demonstrated  that  when  aprotic  RTILs  (having 
only  alkyl  groups  bound  to  the  cation  ring  nitrogens)  are  heated 
in  a  vacuum,  recovery  of  the  vapor  produced  the  same  RTIL 
with  little  or  no  degradation  for  several  families  of  RTILs. 
Vaporization  of  mixtures  of  RTILs  showed  that  one  component 
was  successfully  enriched  upon  fractional  distillation  due  to  the 
differences  in  vapor  pressure  of  one  RTIL  versus  another.7  The 
thermal  stability  of  certain  families  of  RTILs  was  attributed  to 
the  weak  nucleophilicity  (tendency  to  react  with  carbon)  of  the 
anion.  Initially,  the  vaporization  mechanism  was  proposed  to 
be  via  cluster  formation  in  the  gas  phase.5  More  recent  studies 
using  photoionization,17  line  of  sight  mass  spectrometry,18-20 
Lourier  transform  ion  cyclotron  resonance,6,7  and  field  desorp¬ 
tion/ionization21  mass  spectrometric  techniques  suggest  that 
vaporization  of  these  thermally  stable  species  is  via  neutral  ion 
pair  formation  [A+B-]  in  the  gas  phase,  or,  in  the  case  of 
dicationic  ionic  liquids,  via  dication— dianion  [A2+(B-)2]  forma¬ 
tion.20  However,  aprotic  RTILs  with  strongly  nucleophilic  anions 
can  react  when  heated,  and  the  recovery  of  pure  RTILs  by 
vacuum  distillation  may  not  be  feasible. 

Theoretical  modeling  of  RTILs22  has  been  of  great  interest 
in  order  to  predict  such  important  properties  as  heats  of 
formation,  heats  of  vaporization,11,23  melting  and  boiling  points, 
viscosity,  and  thermal  decomposition  mechanisms.24  Dynamics 
calculations  have  indicated  that  vaporization  to  produce  an  intact 
ion  pair  is  energetically  favored  over  ion  cluster  formation 
(multiple  cation— anion  pairs),  where  larger  clusters  have  higher 
heats  of  vaporization.25  The  ionic  nature  of  these  liquids  has 
made  predicting  these  properties  difficult  and  very  few  experi¬ 
mental  data  exist  to  confirm  theoretical  heats  of  formation  and 
heats  of  vaporization  values,  quantities  that  are  essential  for 
evaluating  RTIL  performance  as  propellants. 

The  implication  that  RTILs  evaporate  as  intact  ion  pairs  is 
based  on  detection  of  the  intact  cation  by  mass  spectrometry 
using  photoionization  or  electron-impact  ionization,  as  well  as 
photoelectron  spectroscopy  and  calculation.  It  is  accepted  that 
photoionization  or  electron  impact  ionization  of  the  neutral  ion 
pair,  followed  by  dissociation,  leads  to  a  cation  and  a  neutral, 
and  the  cation  is  then  detected  by  the  mass  spectrometer: 

[A+B  ]  -\-  Jw  [A+B]  +  e  A+  +  B  +  e 

(2) 

Direct  detection  of  the  ion  pair  has  been  elusive,  most  likely 
because  of  the  use  of  relatively  high  energies  of  the  ionizing 
photons  (23.2  eV)  or  electrons  (10—100  eV)  versus  the  energy 
required  to  dissociate  the  ion  pair: 

[BMIM]+[dca]“  —  BMIM+  +  dca' 

A H  =  3.59  eV  at  298  K  (G3MP2)10  (3) 

where  [BMIM]+[dca]-  is  the  RTIL  1 -butyl-3 -methylimidazo- 
lium  dicyanamide.  However,  recent  progress  has  been  made 
toward  direct  detection  of  ion  pairs.21 

In  this  study,  a  first  report  is  provided  for  determinations  of 
the  heats  of  vaporization  of  room  temperature  ionic  liquids  by 
tunable  synchrotron  photoionization  detection.  Several  ionic 
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Figure  2.  Ionic  liquid  effusive  beam  source. 


liquids,  V-butyl-V-methylpyrrolidinium  bistrifluorosulfonylim- 
ide,  [pyrri4]+[NTf2]-,  V-butyl-V-methylpyrrolidinium  dicyana¬ 
mide,  [pyrri4]+[dca]-,  and  1 -butyl-3 -methylimidazolium  dicy¬ 
anamide,  [BMIM]+[dca]-,  are  studied  by  means  of  an  effusive 
vapor  source  coupled  to  a  vacuum  ultraviolet  photoionization 
time-of- flight  mass  spectrometer  (PI-TOLMS).  An  analysis  of 
the  strengths  of  the  photoionization  method  is  made,  which 
include  the  ability  to  study  the  heats  of  vaporization  of  the  ionic 
liquid  molecules  at  the  tunable  thresholds  for  dissociative 
photoionization,  good  comparisons  between  the  A ifvap  results 
from  several  fragment  masses,  and  an  ability  to  assess  the 
temperature  dependence  of  the  photoionization  cross  sections. 
Lurthermore,  the  comparison  of  experimental  photoion  appear¬ 
ance  energies  with  theory  supports  a  dissociative  photoionization 
mechanism.  To  our  knowledge,  this  study  is  the  first  reported 
measurement  of  the  heat  of  vaporization  of  [pyrr14]+[NTf2]-. 
The  heats  of  vaporization  of  [BMIM]+[dca]-  and  [pyrri4]+[dca]- 
obtained  here  are  in  reasonable  agreement  with  previously 
reported  experimental  measurements. 

Experimental  Section 

The  experiments  are  performed  on  the  Chemical  Dynamics 
Beamline,  a  tunable  VUV  beamline  9.0.2. 3  at  the  Advanced 
Light  Source  in  Berkeley,  California.  Product  masses  are 
detected  as  a  function  of  the  vaporization  temperature  of  the 
ionic  liquid  and  tunable  photoionization  wavelength.  The  PI- 
TOLMS  repeller  is  pulsed  at  10  kHz,  and  typically,  data  are 
accumulated  for  200  000  pulses  per  mass  spectrum.  Lor 
[pyrri4]+[dca]-,  200  000  pulses  are  accumulated  at  12.0  eV 
photon  energy  at  473,  483,  and  493  K.  Lor  [pyrri4]+[NTf2]- 
and  [BMIM]+[dca]-,  photoionization  efficiency  (PIE)  data  were 
acquired  as  the  VUV  photon  energy  is  scanned  from  8.0  to  15.0 
eV  in  0.2  eV  steps  (photon  energy  resolution  was  0.025  eV). 
The  effusive  beam  is  oriented  counter  propagating  to  the  VUV 
beam.  In  this  first  study,  PIE  data  are  taken  over  a  limited 
temperature  range,  at  468  and  483  K  for  [BMIM]+[dca]-  and 
at  583  and  613  K  for  [pyrri4]+[NTf2]-.  In  addition,  data  for 
[pyrri4]+[NTf2]-  are  also  accumulated  at  12.6  and  15.0  eV  for 
longer  signal  averaging  of  9  million  pulses  per  spectrum  at  both 
583  and  613  K. 

The  RTIL  source  is  a  modified  version  of  a  source  described 
previously.17  Briefly,  the  source  uses  a  Pyrex  reservoir  that  has 
a  narrow  outlet  tube  (Figure  2).  Each  region  is  surrounded  by 
a  machined  aluminum  block.  Each  block  is  heated  separately 
by  four  cartridge  heaters  embedded  evenly  around  the  cell  axis. 
The  heaters  in  each  region  are  controlled  by  separate  power 
supplies  (0—40  W  per  region),  and  thermal  equilibrium  within 
the  source  is  reached  before  PIE  curves  are  measured.  The 
temperatures  in  each  region  are  measured  using  previously 
calibrated  type  K  thermocouples,  and  from  the  uniformity  of 
the  heated  regions  and  thermocouple  readings,  the  temperature 
values  reported  here  are  assessed  to  have  an  error  of  ±3  K. 
Although  the  temperatures  of  the  reservoir  and  outlet  regions 
can  be  controlled  separately,  for  heat  of  vaporization  determina- 
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Figure  3.  Mass  spectrum  of  [pyrr,4]+[NTf2]-  at  15.0  eV  photon  energy 
and  a  source  temperature  of  613  K. 


Figure  4.  Photoionization  efficiency  (PIE)  curves  of  [pyrri4]+[NTf2] 
for  mass  142  at  583  and  613  K  and  mass  84  at  613  K. 


tions,  the  temperature  of  the  reservoir  region  is  set  to  be  equal 
to  the  outlet  region.  In  order  to  examine  pyrolysis  of  ion  pairs 
in  the  gas  phase,  the  temperature  of  the  outlet  region  can  be  set 
at  a  higher  temperature  than  the  reservoir  region,  and  pyrolysis 
experiments  will  be  described  in  a  future  publication.  The 
experimental  temperatures  are  selected  by  performing  dif¬ 
ferential  scanning  calorimetry  (DSC)  on  the  ionic  liquids  prior 
to  the  photoionization  experiments.  The  temperatures  are 
selected  to  be  well  below  the  decomposition  temperatures  of 
the  ionic  liquids  to  avoid  pyrolysis. 

The  raw  mass  spectra  are  corrected  for  fluctuations  in  the 
power  of  the  light  source  by  monitoring  the  power  at  a  VUV 
photodiode  whose  wavelength  response  curve  had  previously 
been  characterized.  The  synchrotron  is  operated  in  constant  top 
off  mode,  and  the  ring  current  is  kept  constant  at  500  mA.  Low 
source  temperature  spectra  are  scaled  and  subtracted  from  the 
high  temperature  data  to  remove  any  background  signals  not 
due  to  the  RTIL.  Figure  3  shows  a  typical  corrected  mass 
spectrum.  PIE  curves  are  determined  by  integrating  the  peak 
areas  of  a  given  mass  in  each  photon  energy  mass  spectrum. 
Appearance  energies  were  determined  by  a  linear  regression  of 
the  first  five  points  above  the  signal-to-noise  ratio,  with  an 
estimated  uncertainty  of  ±0.2  eV. 

1 -Butyl-3 -methylimidazolium  dicyanamide  (>98%)  and 
V-butyl-V-methylpyrrolidinium  dicyanamide  (>99%)  are 
obtained  from  EMD  Chemicals/Merck.  V-Butyl-V-meth- 
ylpyrrolidinium  bistrifluorosulfonylimide  (>98%)  is  pur¬ 
chased  from  Aldrich.  In  order  to  remove  volatile  impurities 
present  in  the  sample,  the  RTIL  samples  are  heated  in  high 
vacuum  (base  pressure  ~1CT6  Pa)  for  12  h  prior  to  the 
experiments  as  follows:  [pyrr14]+[NTf2]_  at  453K, 
[pyrri4]+[dca]-  at  378  K,  and  [BMIM]+[dca]“  at  358  K. 

Ab  initio  calculations  are  performed  using  Gaussian  03  W26 
at  the  HF/6-31±G(d,p)  or  B3LYP/6-31±G(d,p)  level  of  theory 
to  preoptimize  geometries.  Final  optimization  and  energy 
calculations  are  performed  at  the  MP2/6-31±G(d,p)  and 
M06/27  6-31±G(d,p)  levels  of  theory  using  GAMESS. 28,29 
Resulting  MP2  and  M06  energies  at  0  K  are  corrected  for 
unsealed  zero-point  vibrational  energies.  On  the  basis  of 
previous  calculations  of  RTIL  ionization  processes  using  HF/ 
6-31±G(d,p),  B3LYP/6-3 1  ±G(d,p),  and  MP2/6-31±G(d,p) 
levels  of  theory,  MP2/6-31±G(d,p)  is  selected  as  the  preferred 
method  due  to  good  agreement  with  experimental  ionization 
potentials  and  ionization  trends.  The  MP2/6-31±G(d,p)  level 


of  theory  is  assigned  an  uncertainty  of  ±0.1  eV  based  on  our 
previous  calculations.  Due  to  complications  with  MP2  for 
[pyrri4]+[dca]-  explained  below,  M06  was  used  to  calculate  the 
ionization  energies  for  [pyrri4]+[dca]“. 

Results  and  Discussion 

As  noted,  a  mass  spectrum  of  [pyrri4]+[NTf2]_  taken  at  613 
K  and  15.0  eV  photon  energy  is  shown  in  Figure  3,  and  the 
PIE  curves  for  mass  142  at  583  and  613  K  and  mass  84  at  613 
K  are  shown  in  Figure  4.  Mass  142  represents  the  [pyrri4]+ 
parent  cation  of  the  ionic  liquid  ion  pair  in  the  gas  phase  (eq 
2).  The  appearance  of  mass  84  is  due  to  fragmentation  of  the 
internally  hot  [pyrri4]+  cation,  which  can  be  attributed  to  C4Hi0 
loss.  From  these  curves,  the  appearance  energies  of  the  resulting 
intact  cation  and  fragment  cations  can  be  determined,  allowing 
insight  into  the  energetics  of  ionization  of  the  ion  pair  and 
possible  ion  fragmentation  pathways.  The  observed  appearance 
energies  for  masses  142  and  84  are  9.4  ±  0.2  and  10.0  ±  0.2 
eV,  respectively,  which  are  in  excellent  agreement  with  their 
corresponding  MP2  appearance  energies  of  9.2  ±0.1  and  10.1 
±0.1  eV.  The  thermal  contribution  due  to  the  small  temperature 
differences  in  these  experiments  (<30  K)  should  only  amount 
to  a  few  millielectron  volts  at  the  temperatures  involved,  causing 
a  shift  to  lower  appearance  energy  due  to  hot  bands  that  is 
smaller  than  the  resolution  of  our  experiment. 

By  comparing  the  intensities  of  mass  142  peaks,  which  are 
proportional  to  the  number  densities  of  the  ion  pair  in  the  RTIL 
vapor,18  at  two  known  source  temperatures  T\  and  T2  and 
applying  the  Clausius— Clapeyron  equation,  the  heat  of  vapor¬ 
ization  in  this  temperature  range  is  directly  determined: 

A Hmp  =  -R  ln(/2//1)[(7’17'2)/(7'2  -  7,)]  (4) 

where  /  is  the  measured  mass  142  peak  intensity  and  R  is  the 
gas  constant.  Over  the  temperature  range  of  this  experiment, 
the  photoionization  cross  section  at  a  given  photon  energy  is 
considered  to  be  constant.  The  branching  fractions  from  the 
ionized  ion  pair  to  intact  cation  or  cation  fragments  are 
considered  to  be  constant  as  well. 

The  heat  of  vaporization  of  [pyrri4]+[NTf2]_  is  determined 
using  eq  4  and  the  data  in  the  following  ways:  (1)  using  the 
sum  of  the  intensities  of  all  the  points  along  the  mass  142  tunable 
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TABLE  1:  Heats  of  Vaporization  Data" 

2010 

Chambreau  et  al. 

[pyrr14]+[NTf2]- 

Atfvap(PIE) 

A//vap(12.6  eV) 

A//vap(15.0  eV) 

A//Vap(average) 

mass  142 

180.3 

156.6 

174.6 

170.5 

mass  84 

167.9 

155.3 

165.4 

162.9 

average  of  six  values  = 

166.7 

correction  to  298.15  K  = 

28.2 

A//vap(298.15)  = 

195  ±  19 

[pyrri4]+[dca]~ 

A//vap(483/473  K) 

A//vap(493/483  K) 

A//vap(493/473  K) 

A//Vap(average)  AHvap(slope) 

mass  84,  12.0  eV 

145.2 

157.1 

151.1 

151.1  151.0 

correction  to  298.15  K  = 

19.5 

A//vap(298.15)  = 

171  ±  12 

[BMIM]+[dca]  A//vap(mass  82)  A//vap(mass  97)  A//vap(mass  124)  A//vap(mass  137)  A//Vap(average) 

PIE  483/468  K  148.4  157.4  162.4  154.8  155.8 

correction  to  298.15  K  =  18.7 

A//vap(298.15)  =  174  ±12 

a  All  values  are  in  kJ/mol.  Uncertainties  in  the  present  work  are  twice  the  standard  deviation. 


TABLE  2:  Experimental  and  Literature  Values  for  Heats  of  Vaporization  and  Temperatures  of  Decomposition 

RTIL 

T  (K)  range 

this  work  A//vap 
(kJ  mol-1)  (average) 

ACPAT  (kJ  mol-1) 

this  work  A//vap 
(kJ  mol-1)  (298.15  K) 

literature  Affvap 
(kJ  moP1)  (298.15  K) 

Tdccom position  (K) 

[BMIM]+[dca]“ 

468-483 

156  ±  12 

18.7 

174  ±  12 

157.2  ±  1.1* 

513* 

[pyrri4]+[dca]_ 

473-493 

151  ±  18 

19.5 

171  ±  12 

161.0  ±2.0C 

535“ 

[pyrr14]+[NTf2]- 

583-613 

167  ±  19 

28.2 

195  ±  19 

612“ 

a  Reference  10.  ^Reference  31.  c  Reference  18.  ^References  16  and  31. 


wavelength  PIE  curve  (36  points  with  7.2  million  pulses  total 
per  curve)  at  T\  and  T2 ,  (2)  by  individually  comparing  the 
intensities  at  T\  and  T2  of  mass  142  at  one  specific  wavelength, 
12.6  eV,  and  (3)  by  analyzing  the  results  also  at  15.0  eV.  A 
longer  signal  averaging  condition  of  9  million  pulses  is  used 
for  the  last  two  cases.  T\  and  T2  are  583  and  613  K,  respectively. 
A  similar  analysis  is  performed  for  mass  84  also.  Table  1 
summarizes  the  experimental  A 77vap  values  obtained.  The 
experimental  A77vap  values  have  been  corrected  to  298.15  K 
using  a  A Cp  value  of  —105.4  J#K-1* mol-1  for  [BMIM]+[dca]“ 
and  [pyrr14]+[dca]~  18  and  a  value  of  —94  J*K-1 ‘mol-1  for 
[pyrr14]+[NTf2]-.10 

A//vap(298.15  K)  =  Ai7vap(7)  +  ACpAT  (5) 

Uncertainties  in  the  A Cp  values  were  not  reported  in  ref  10, 
but  their  influence  on  the  temperature  corrected  A 7/vap  values 
was  estimated  to  be  negligible  within  experimental  uncertainties, 
and  therefore  they  are  not  included  here.  The  uncertainties  in 
Table  1  represent  twice  the  standard  deviation  of  the  measure¬ 
ments.  The  larger  uncertainties  in  these  measurements  are  the 
result  of  having  fewer  data  points  across  similar  temperature 
ranges  for  determining  the  heat  of  vaporization  versus  the  more 
extensive  maps  versus  temperature  in  the  referenced  works. 

Good  agreement  for  A77vap  calculated  using  masses  142  and 
84  (170.5  vs  162.9  kJ/mol,  Table  1)  indicates  the  possibility  of 
calculating  A 7/vap  from  mass  fragments  of  the  cation,  and  mass 
84  is  used  successfully  to  determine  the  A 7/vap  of  [pyrri4]+[dca]~. 
One  possibility  why  mass  84  is  usable  for  determining  the  heat 
of  vaporization  for  the  [pyrrol +-containing  RTILs  may  be  due 
to  its  increased  stability  from  possible  resonance  structures  of 
the  cation.18  This  stability  also  could  explain  the  significant 
fragmentation  of  the  intact  cation  to  mass  84  under  such  soft 
ionization  conditions.  An  additional  Clausius— Clapeyron  cal¬ 
culation  was  performed  by  plotting  In  I  versus  1  IT  of  mass  84 
at  three  temperatures,  namely,  473,  483,  and  493  K.  The  heat 


of  vaporization  determined  by  the  slope  of  the  straight  line 
through  these  points  agrees  well  with  the  average  value 
(A77vap(  slope)  vs  A77vap(average),  Table  1).  Mass  84  analysis 
for  the  heat  of  vaporization  of  [pyrrol +[dca]“  agrees  reasonably 
well  with  the  reported  literature  value  (Table  2).  For 
[BMIM]+[dca]~,  masses  82,  97,  124,  and  137  are  used  to 
calculate  the  heat  of  vaporization,  and  the  reported  heat  of 
vaporization  value  is  slightly  higher  than  the  literature  value 
(Table  1).  The  successful  determination  of  A//vap  from  fragment 
masses  is  important  because  for  some  RTILs  the  intact  cation 
may  not  be  detected.  In  fact,  only  a  very  weak  parent  cation 
ion  peak  at  mass  142  is  detected  for  [pyrri4]+[dca]~  and  no 
parent  peak  at  mass  139  is  detected  for  [BMIM]+[dca]“.  The 
present  A 77vap  of  [pyrr14]+[NTf2]“  is  consistent  with  heats  of 
vaporization  determined  for  other  NTf2“  based  RTILs.6’12’13’23’25’30 
The  temperatures  used  in  these  experiments  are  well  below  the 
reported  onset  decomposition  temperatures  (nonzero  DSC  slope) 
for  these  species  (Table  2). 

To  assess  the  temperature  dependence  of  the  photoionization 
cross  section,  the  relative  shapes  of  the  PIE  curves  obtained  at 
the  experimental  reservoir  temperatures  are  compared.  Figure 
5a  shows  the  comparison  of  the  [pyrri4]+[NTf2]“  mass  142  PIE 
curves  at  583  and  613  K.  In  order  to  try  to  force  the  curve  shapes 
to  match,  rather  than  simply  scaling  the  PIE  curve  at  583  K  to 
the  PIE  curve  at  613  K  so  that  the  highest  points  match,  the 
PIE  curve  at  583  K  is  scaled  iteratively  to  fit  to  the  613  K  data 
by  a  least- squares  fit  where  the  sum  of  the  squares  of  the 
difference  between  two  points  at  a  fixed  photon  energy  is 
minimized: 

(d/dS')2(S7583  -  /613)2  =  0  (6) 

By  this  method,  if  the  scaled  PIE  curves  are  identical,  then 
2(S7583  —  /613)2  would  necessarily  equal  zero,  and  S  =  (/6 13/ 
/5 83).  However,  if  the  curves  do  not  match  exactly,  then  2(5/58 3 
—  76i3)2  >  0.  In  Figure  5a  and  b,  for  [pyrri4]+[NTf2]-,  the  high 
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Figure  5.  Least-squares  fit  analyses  of  the  PIE  curves  for  (a)  [pyrri4]+[NTf2]  of  mass  142  at  583  K  relative  to  613  K,  (b)  [pyrrol +[NTf2]  of  mass 
84  at  583  K  relative  to  613  K,  (c)  [BMIM]+[dca]~  of  mass  97  at  468  K  relative  to  483  K,  and  (d)  [BMIM]+[dca]~  of  mass  137  at  468  K  relative 
to  483  K. 


temperature  PIE  curves  for  both  mass  142  and  mass  84  appear 
to  have  shifted  slightly  to  higher  intensities  between  10  and 
13.6  eV  ( 57583  <  /6 i3)  and  to  lower  intensities  above  13.6  eV 
(SI >  hn)  relative  to  the  low  temperature  PIE  curves.  Due  to 
this  slight  shift  in  the  PIE  shape  at  higher  temperature,  in  the 
region  where  the  PIE  intensity  increases  with  temperature,  *S7583 
<  /6i3  (10—13.6  eV),  the  Clausius— Clapeyron  approach  sys¬ 
tematically  overestimates  the  heat  of  vaporization.  In  the  region 
where  the  PIE  intensity  decreases  with  temperature  above  13.6 
eV  where  Sis 83  >  hn,  the  Clausius— Clapeyron  calculation  will 
underestimate  the  heat  of  vaporization.  Because  there  are  more 
points  in  the  PIE  summed  calculation  where  5758 3  <  h\3  (8—13.6 
eV,  29  points)  than  where  S7583  >  I613  (13.8—15  eV,  7  points), 
the  PIE  summed  result  from  the  Clausius— Clapeyron  calculation 
will  systematically  overestimate  the  actual  heat  of  vaporization, 
although  these  effects  will  tend  to  cancel  somewhat.  If,  instead 
of  determining  the  heat  of  vaporization  of  [pyrr  i4]+[NTf2]  using 
the  ratio  of  summed  PIE  intensities  (hn/Isss  —  6.17,  A Hwap  = 
180.3  kJ/mol,  Table  1)  of  mass  142  in  the  Clausius— Clapeyron 
equation  one  uses  the  curve  fitting  scaling  factor  S  =  5.29  from 
Figure  5a,  the  result  is  A i/vap  =  165.0  kJ/mol,  a  reduction  in 
the  calculated  value  of  A//vap  by  8.9%.  A  similar  treatment  for 
mass  84  (Figure  5b)  and  for  masses  97  and  137  of 
[BMIM]+[dca]“  in  Figure  5c  and  d  results  in  overestimates  of 
1.1,  6.3,  and  1.1%,  respectively.  The  similar  shifts  in  the  PIE 
curves  in  Figure  5  indicate  this  is  a  systematic  phenomenon. 
On  the  basis  of  A 7/vap  calculated  from  the  measured  PIE  curve 


sum  values  and  comparison  with  the  scaling  factor  analysis 
which  takes  into  account  the  slight  temperature  shift  in  PIE 
curves,  the  A Hvav  values  determined  by  PIE  summed  data  likely 
contain  a  systematic  overestimation  of  1  —  10%,  depending  on 
how  much  the  PIE  curve  shifts  at  the  higher  temperature. 
Although  the  A 7/vap  values  determined  by  PIE  were  included 
in  the  reported  average  value  for  [pyrri4]+[NTf2]“  in  Table  1, 
the  uncertainty  (±19  kJ/mol,  twice  the  standard  deviation)  is 
9.7%  of  this  value  and  encompasses  this  systematic  error.  An 
improvement  in  the  accuracy  of  this  method  will  be  to  select  a 
photon  energy  where  the  scaled  PIE  curves  cross  at  S  =  /2//i 
and  competition  between  photoionization  cross  section  and  ion 
fragmentation  cancels  each  other  out.  For  example,  monitoring 
the  vapor  density  as  a  function  of  temperature  at  ~13.6  eV  and 
determining  A 77vap  from  the  slope  of  a  plot  of  In  I  versus  1/T18 
would  increase  both  the  accuracy  and  precision  of  this  method. 

The  PIE  curve  shift  at  higher  temperature  could  be  due  to 
several  factors.  Competition  between  increased  photoionization 
cross  section  and  increased  ion  fragmentation  at  high  temper¬ 
atures  may  account  for  this  observed  shift.  Below  13.6  eV,  a 
larger  photoionization  cross  section  at  high  temperatures  could 
explain  the  shift  to  higher  intensities,  and  above  13.6  eV, 
increased  ion  fragmentation  reduces  the  PIE  intensity.  A  direct 
comparison  of  the  branching  ratios  of  masses  84/142  at  two 
different  temperatures  is  shown  in  Figure  6.  Above  ~12  eV, 
the  increase  in  branching  ratio  at  higher  temperature  indicates 
increased  ion  fragmentation  with  increasing  temperature. 
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Figure  6.  Mass  84/142  branching  ratios  for  [pyrri4]+[NTf2]  at  583 
and  613  K  as  a  function  of  photon  energy. 

There  are  several  advantages  for  using  photoionization  versus 
electron  impact  ionization  for  measuring  heats  of  vaporization 
for  RTILs.  The  narrow  energy  spread  of  the  light  source  (0.025 
eV)  and  soft  ionization  conditions  minimize  the  effect  of 
fragmentation  of  the  parent  cation  and  allow  for  appearance 
energies  to  be  determined  accurately.  Typical  electron  impact 
ionization  energies  are  around  70  eV,  where  the  electron  impact 
ionization  cross  section  curve  is  large  and  the  shape  of  the  curve 
is  flat  for  many  atomic  and  molecular  species,  but  molecular 
fragmentation  to  smaller  species  can  be  extensive.  In  this  work, 
photoionization  measurements  yield  a  branching  ratio  of  masses 
84/142  of  nearly  1  for  [pyrri4]+[NTf2]_,  indicating  the  small 
extent  of  fragmentation  of  the  [pyrr]+  under  the  present  soft 
ionization  conditions  versus  electron  impact  ionization  of 
[pyrri4]+[dca]“  with  a  branching  ratio  of  masses  84/142  that  is 
about  an  order  of  magnitude  higher.18  Although  fragmentation 
of  the  parent  ion  can  be  minimized  by  lowering  the  electron 
impact  energy,  the  electron  impact  ionization  cross  section 
changes  rapidly  as  a  function  of  electron  energy  near  the 
ionization  threshold,31,32  and  may  be  highly  susceptible  to 
thermal  effects  near  the  ionization  threshold.  Also,  the  energy 
spread  of  a  typical  electron  gun  is  limited  to  >0.5  eV,  limiting 
the  energy  resolution  in  threshold  measurements.  Therefore, 
determining  appearance  energies  of  ions  by  electron  impact 
ionization  is  difficult  and  has  large  associated  uncertainties. 

From  the  accurate  determination  of  photoion  appearance 
energies  by  the  present  method,  comparison  with  theoretical 
energy  calculations  can  be  made.  Ab  initio  calculations  of  the 
energetics  involved  in  the  dissociative  photoionization  processes 
(described  below)  are  in  excellent  agreement  with  the  experi¬ 


mentally  determined  photoion  appearance  energies.  The  mea¬ 
sured  mass  142  appearance  energy  in  Figure  4  (9.4  d=  0.2  eV) 
agrees  well  with  the  ab  initio  value  of  9.2  ±0.1  eV  calculated 
at  the  MP2/6-31±G(d,p)  level  of  theory  (0  K,  ZPVE  corrected) 
for  the  [pyrri4]+[NTf2]“  ion  pair  (eq  7).  An  appearance  energy 
of  10.1  ±  0.1  eV  (experimental  value  =  10.0  ±  0.2  eV)  was 
calculated  for  mass  84,  assuming  elimination  of  C4Hi0  from 
the  pyrri4+  cation  to  form  the  resonance  stabilized  cationic 
structure  postulated  in  a  previous  publication  (eq  8). 18  The 
electron  affinity  for  the  NTf2*  radical  is  calculated  to  be  5.7  eV 
at  the  MP2/6-31±G(d,p)  level  of  theory. 
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Calculation  of  the  adiabatic  ionization  potentials  of  dicyanamide 
and  [pyrri4]+[dca]“  at  the  MP2  level  was  precluded  due  to  the 
inadequacy  of  a  single  configuration  reference  wave  function 
for  the  neutral  dicyanamine  species.  (The  need  for  a  multicon- 
figurational  reference  in  dicyanamine  is  evident  from  (a)  the 
high  degree  of  spin  contamination  present  in  the  unrestricted 
Hartree— Fock  (UHF)  wave  function  (( S 2>  =  1.303,  compared 
to  0.5  for  a  pure  spin  doublet)  and  (b)  the  presence  of  multiple 
unphysical  natural  orbital  occupation  numbers  (see  Gordon, 
M.  S.;  Schmidt,  M.  W.;  Chaban,  G.  M.;  Glaesemann,  K.  R.; 
Stevens,  W.  J.;  Gonzalez,  C.  J.  Chem.  Phys.  1999,  110 ,  4199) 
obtained  from  the  second-order  Z-averaged  perturbation  theory 
(ZAPT(2))  relaxed  density  matrix.)  However,  the  electron 
affinity  for  the  dicyanamide  radical  has  been  determined 
experimentally  to  be  4.135  eV.33  Using  this  value  along  with 
the  theoretical  values  in  Table  3,  the  appearance  energy  for  mass 
84  in  eq  10  is  estimated  to  be  8.73  ±0.1  eV,  and  is  in  excellent 


TABLE  3:  Ab  initio  Energies  (in  Hartrees)  of  Proposed  Species  Involved  in  Dissociative  Photoionization  of  [pyrri4]+[NTf2]  and 
[pyrri4]+[dca] 


species 

E0 

zero  point  energy 

total  H0 

[pyrri4]+[NTf2]~  ion  pair 

-2232.073761 

0.348228 

-2231.725533 

[pyrri4]+[NTf2#]  ion  pair  cation 

-2231.758743 

0.346514 

-2231.412229 

[pyrri4]+  cation 

-408.215636 

0.293093 

-407.922543 

[NTf2]~  anion 

-1823.729958 

0.053730 

-1823.676228 

NTf2*  neutral  radical 

-1823.518844 

0.052807 

-1823.466037 

[pyrri4]+[dca]~  ion  pair 

-648.205863 

0.314189 

-647.891674 

[pyrri4]+[dca#]  ion  pair  cation  (MP2  vertical) 

-647.924455 

0.314189 

-647.610266 

[pyrri4]+[dca]~  ion  pair  (M06) 

-649.732310 

0.306353 

-649.425957 

[pyrri4]+[dca#]  ion  pair  cation  (M06  vertical) 

-649.444280 

0.306353 

-649.137927 

[pyrri4]+[dca#]  ion  pair  cation  (M06  adiabatic) 

-649.448255 

0.034323 

-649.143932 

[dca]~  anion 

-239.854100 

0.020016 

-239.834084 

mass  84  cation 

-250.263014 

0.151086 

-250.111928 

butane 

-157.912441 

0.135673 

-157.776768 
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agreement  with  I  he  experimental  appearance  energy  for  maws 
84  in  the  [pyrrol +[  do  a]-  dissociative  photoionization  of  8,8  _L 
0,2  eV,  The  estimated  appearance  energy  of  The  intact  [pyiTi4]+ 
cation  in  eq  9  is  7.81  eV.  However,  the  experimental  signal- 
tO-noi$e  ratio  for  this  species  was  not  sufficient  to  determine 
its  appearance  energy,  although  ir  was  detectable  at  higher 
energies-.  From  M06  calculations*  the  vertical  and  adiabatic 
ionization  potentials  of  [pyrrt4j+|dcaj“  are  7,8  and  7,7  eV, 
respectively,  and  the  M06  vertical  ionization  potential  agrees 
well  with  the  MP2  vertical  ionization  potential  of  7.7  eV, 
( )  veral  1 ,  I  he  agrees n  en  t  hel w  een  1 1 1  eory  an  t  i  exp  er  n n  edit  s  upport  s 
the  dissociative  photoionization  of  isolated  cation— anion  pairs 
in  the  gas  phase. 

Conclusions 

Heats  of  vaporization  for  several  room  temperature  ionic 
liquids  have  been  determined  by  effusive  beam  VUV  photo¬ 
ionization  mass  spectrometry.  The  photoionization  cross  section 
is  considered  to  be  temperature  independent  when  using  die 
PIE  curves  for  the  Clausius— Clap eyron  analysis.  We  estimate 
a  positive  systematic  error  of  up  to  10%  for  this  assumption. 
After  thermal  correction  to  298.15  K,  the  A values  of 
[BMIM]+[dca]“  and  [pyrrol"1" [deal-  reported  here  arc  in  reason¬ 
able  agreement  with  l he  previous ty  reported  experimental  values. 
For  [pynv.]+[NTf£T*  A//V^<298.15  K)  is  calculated  to  be  195 
d:  19  k,T  mol-1,  which  has  not  been  reported  previously. 
Detection  of  the  intact  [pyrrol*  cation  is  indicative  of  the 
presence  of  neutral  cation  anion  pairs  in  the  effusive  beam 
which  dissociate  to  form  [pynrJ*  and  neutral  NTfj-  at  low  VUV 
photon  energy.  A  second  dissociation  pathway  detected  is  the 
efficient  fragmentation  of  [pynu] '  to  produce  mass  84, 
Together,  these  experimental  heals  of  vaporization,  photoion¬ 
ization  appearance  energies,  and  ab  initio  calculations  support 
a  vaporization  mechanism  for  RTILs  as  intact  cation— anion 
pahs. 
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Abstract 

lliis  paper  deals  with  the  gas-phase  thermodynamic  properties 
of  endothermic  compounds  potentially  formed  during  mouame- 
thylhydrazinc  ( M MH  V ni ti ogc n  tetroxide  (N  lO)  hypcrgolic  re¬ 
act  iv  it  y.  The  standard  enthalpies  ul  format  ion  at  298.15  K  are 
detemiiiied  by  means  of  quantum  chemistry  calculations  along 
with  protocols  developed  for  these  compounds.  The  resultant 
data,  currently  previously  unavailable  for  almost  all  of  these 
compounds,  are  potentially  critical  to  the  modeling  of  combustion 
chemistry  of  this  bipropcllant  combi  nation. 

Keywords:  Standard  Enthalpies  of  Formation,  Monomethyl- 
hydrazinc.  Nitrogen  Ictroxidc,  U  asymmetrical  Dimethyl 
Hydrazine 

1  Introduction 

Nitrogen  tetroxide  (NTO)/monomethylhydraztno 
(MMH),  NTO/u  ns  ym  metrical  dimethyl  hydrazine 
(UDMII)  and  NTO/Aerozine  50  (MMi  I/UDMH  mixture) 
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are  chemical  systems  used  for  in-space  liquid  propulsion. 
These  systems  are  very  interesting  because  they  are  able  to 
auto  ignite  at  low  temperature  (room  temperature  and 
below),  he.  lliey  ignite  without  any  external  ignilion  device. 
This  phenomenon  is  called  hypergolicily.  From  a  funda¬ 
mental  research  point  of  view,  much  more  work  is  needed  to 
understand  all  the  features  of  hypcrgolic  ignition.  Fuels  are 
stored  in  the  liquid  state  but,  as  usual  for  liquids,  ignition 
(hypcrgolic  or  not)  occurs  in  the  gas-phase  above  the  liquids. 
This  has  been  recently  shown  for  a  hypcrgolic  system  by 
Alfano  et  al.  [I|.  Gas-phase  species  and  eonde used-phase 
species  form  during  the  reaction  between  NTO  and  MMH 
or  UDMH*  Experimental  observations  show  that  con¬ 
densed  species  form  when  the  reaction  does  not  lead  to 
ignition  and  that  condensed  species  form  when  hypcrgolic 
ignition  occurs,  i.e.  during  the  ignition  delay  |2j.  The 
chemical  nature  of  these  compounds  is  not  clearly  estab¬ 
lished,  despite  numerous  studies.  In  fact,  the  NTO/MMH 
chemistry  depends  on  the  experimental  conditions:  equiv¬ 
alence  ratios,  temperature,  and  pressure.  liquid-liquid 
reactions  at  low  temperature  lead  to  condensed  methylky- 
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drazinium  nitrate,  It  can  also  be  the  case  that  methylhy- 
drazinium  azide  forms.  However,  even  if  these  species  can 
play  a  role  in  other  related  topics,  such  as  hard  starts,  their 
role  in  hypergolic  ignition  is  not  completely  clear.  A  detailed 
chemical  kinetic  model  was  proposed  to  interpret  the  gas- 
phase  NTO/MMH  reactivity  1 3].  This  mechanism  includes  a 
sub- mechanism  devoted  to  MMH  thermal  decomposition, 
validated  with  MMH  profiles  observed  in  shock  lube 
experiments  [4]  and  a  sub-mechanism  devoted  to  the 
oxidation  of  MMH  validated  with  MMH/fh  [5]  and 
MMH/O2/H2  [6]  ignition  delays.  The  formation  of  mcthyl- 
diazene  CH3N— NH  and  its  subsequent  thermal  decompo¬ 
sition  together  with  the  exothermic  formation  of  nitrite  and 
nitro  compounds  resulting  from  ter  molecular  recombina¬ 
tions  between  mrthyihydrazyl  radicals  andNCK  explain  gas- 
phase  hypergolic  ignition  and  gas-phase  non-ignition.  It 
remains  uncertain  whether  this  model  is  sufficiently  de¬ 
tailed,  but  it  allows  the  simulation  of  gas-phase  MM1 1/NTG 
self-ignition  at  low  temperature  (ambient  and  below),  for 
any  pressure,  and  this  is  certainly  a  major  achievement.  It  is 
commonly  the  ease  that  only  about  10%  of  the  reactions 
given  in  a  detailed  chemical  kinetic  model  play  a  role  on  l he 
ignition/comb  ust  ion  process  and  l  here  Core  the  writing  of  all 
the  likely  reactions  is  not  always  appropriate.  Sun  and  I  .aw 
[7j  recently  proposed  rate  constants  for  some  elementary 
reactions  needed  to  simulate  the  thermal  decomposition  of 
MMH.  These  reactions  are  not  expected  to  play  a  role  in 
MMH/NTO  hypergolic  ignition,  but  this  remains  to  be 
demonstrated.  Nonnenberg  et  al.  [8J.  and  Prank  etal.  [9] 
used  molecular  dynamics  to  simulate  NTO/MMH  reactivity 
and  found  that  methyldiazene  can  be  assumed  to  be  the 
main  product  during  NTO/MMH  hypergolic  ignition. 
According  to  this  paper,  the  recombination  of  methylhy- 
drazyl  radicals  can  lead  to  di  methyl  tetrazanes,  Indeed  these 
termoleeulai  recombination  reactions  ar  e  exothermic  and 
must  be  considered  in  the  NTO/MMH  delailed  chemical 
kinetic  model.  Thermodynamic  data  needed  for  the  writing 
of  MMH-based  detailed  kinetic  models  are  presented  in 
Ref.  110]  and  [11].  The  thermochemistry'  presented  by  Sun 
and  Laur  ]7]  for  m  ethyl  hydra  zyl  radicals  CH^N'NHj,  CHi 
NHMPand  CH2'NHNH3  is  consistent  with  the  data  pre¬ 
sented  by  Caloire  and  Swiharl  [10).  However,  the  thermo¬ 
chemistry  of  the  dimethyl  tetrazanes  is  not  available  and 
remains  to  be  studied.  McQuaid  et  al.  [12,  13]  also  present 
considerations  for  MMH/lRi  NA  hypergolic  ignition,  which 
partially  validate  the  NTO/MMH  detailed  chemical  kinetic 
model  proposed.  Finally,  the  formation  of  a  number  of 
nitroso,  nitrite,  nitro,  and  nitrate  compounds,  not  previously 
taken  into  account,  must  be  considered  for  completeness, 
Nitroso,  nitro,  nitrile,  and  nitrate  compounds  considered 
here  have  an  —NO  group,  -NG2  group,  — ONO  group,  or 
-0N02  group  attached  to  nitrogen  atom,  respectively. 
Experimental  data  for  such  derivatives  of  MMH  and 
UDMH  are  scarce.  The  thermochemistry  of  these  com¬ 
pounds  is  needed  to  allow  the  construction  of  more  detailed 
chemical  kinetic  models  able  to  simulate  hypergolic  igni¬ 
tion.  The  aim  of  this  work  is  to  present  gas-phase  standard 
enthalpy  of  formation  of  nitroso,  nitro,  nitrite,  and  nitrate 


derivatives  of  MMH  and  1 3DMH,  not  considered  previous¬ 
ly,  and  gas-phase  standard  enthalpy  of  formation  of 
dimethyl l etrazane  compounds.  These  data  are  obtained  by 
using  Iwo  ab  initio  methods,  specifically  derived  and 
validated  for  this  purpose,  one  for  compounds  containing 
NOj,  ONO,  ON  On  groups  and  one  for  general  organic 
compounds. 

2  Computational  Details 

The  gas-phase  standard  enthalpy  of  formation  of  mole¬ 
cule  j  at  298.15  K  can  be  determined  from  the  following 
equation: 

4^..^  =  2625.5 

x  +  ZPEj  l  thermal  corrections  +  V  J 

where  «r  is  the  number  of  atoms  i  in  molecule  j  and  c]  the 
atomic  correction  for  atom  L  Ej  and  Z PE;  denote,  respec¬ 
tively,  the  absolute  electronic  energy  and  zero-point  energy, 
calculated  using  the  Gaussian  9RW  [14]  and  Gaussian  03  [  J  5] 
softwares,  't  he  u  ni is  are  I  la  rt  ree  mo  lecule  1  for  Ej ,  ZPEy  and 
thermal  corrections,  and  Haiti  ee  atom  1  for  c*lt  whereas 
is  in  kJ  mol  J.  Two  classes  of  compounds  are 
considered  in  this  study:  C/H/N /Q  compounds  and  C/H/N 
compounds.  The  method  established  and  validated  for 
nitroso,  nitro,  nitrite,  and  nitrate  compounds,  i.e.  C/H/N/O 
compounds,  is  described  in  Osmont  et  al.  [16].  the  method 
established  and  validated  for  tetrazanes,  l.e.  C/H/N  com¬ 
pounds.  is  described  in  Osmont  et  al.  [17,  18]. 

3  Results 

Experimental  difficulties  in  isolating  these  compounds  for 
calorimetric  Studies  arc  substantial,  and  therefore  estima¬ 
tion  methods  have  to  be  used.  Such  methods  can  be  more  or 
less  sophisticated.  In  this  study,  we  used  two  ab  initio 
methods  with  empirical  atomic  corrections.  The  details  of 
these  methods  are  beyond  the  scope  of  this  paper.  One  of 
these  methods  was  specifically  established  and  validated 
with  experimental  data  known  011  nitro,  nitrite,  and  nitrate 
compounds.  This  method  was  derived  wrilh  data  from  28 
compounds  and  validated  with  16  compounds.  The  average 
absolute  deviation  between  experiments  and  calculations  is 
7.1  kJ  mol  1  for  the  derivation  set  and  10.0  kT  mol  1  for  the 
validation  set.  Table  1  gives  the  comparison  between 
calculated  and  experimental  gas-phase  standard  enthalpies 
of  formation  at  298.15  K  for  some  compounds  included  in 
the  derivation  and  validation  datasets.  The  agreement 
between  experiments  and  calculations  is  good.  Tire  average 
absolute  deviation  for  this  method,  in  general,  is  expected  to 
be  =L  13  kJ  mol  K  In  table  1,  the  average  absolute  deviation 
is  5.4  k.I  mol  1  for  the  nitrite  compounds,  R.ft  kJ  mol  1  for  the 
nitro  compounds  and  11.3  kJ  mol  1  for  the  nitrate  com- 
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Tabic  1.  Comparison  between  experimental  and  calculated  (within  parenthesis)  gas-phase  standard  enthalpy  of  formation  at  298.15  K 
for  nitrite,  nitro,  and  nitrate  compounds  from  the  validation  database.  Units  are  kJ  mol-1. 


-X 

-ONO 

-no2 

-ONO? 

CHjX 

—  65.3  (-65.3) 

-80.7  (-82.0) 

-122.2  (-129.3) 

QHjX 

-108.4  (-102.9) 

-102.5  (-117.1) 

-154.8  (-166.5) 

QH;X 

-118.8  (-128.9) 

-  124.7  (- 138.5) 

-174.0  (-187.9) 

CH,-CHX— CH? 

-133.5  (-141.0) 

-140.2  (-150.2) 

-190.8  (-202.5) 

(CHj)3C— X 

-171.5  (-167.4) 

-177.0  (-180.7) 

No  experimental  data 

Table  2.  Calculated  gas  phase  standard  enthalpies  at  298.15  K  for 
nitToso,  nitrite,  nitro,  and  nitrate  derivatives  of  MMII  and 
UDMHL  Units  are  ki  mol-1. 


-X 

-NO 

-ONO 

-no2 

-ono2 

ch3-nh-nh-x 

197.9 

173.6 

102.1 

142.3 

ch-nx-nh2 

158.2 

173.6 

95.8 

113.8 

ch3-nh-nx. 

342.7 

292.0 

183.8 

172.4 

ch3-nx-nh-x 

278.7 

284.9 

143.9 

167.4 

CH -NX~NX? 

427.6 

387.0 

247.7 

232.6 

(CHjV-N-NH-X 

174.9 

187.9 

103.8 

123.8 

(ch3)2-n-nx2 

319.2 

288.3 

182.0 

177.0 

Table 3.  Gas  phase  standard  enthalpy  of  formation  at  298  K  of 
nitro  and  nitrite  derivatives  of  MMII  calculated  by  McQuaid  et  al. 
[12. 13].  Units  arc  U  mol  \ 


-X 

-ONO 

-  no2 

CHj-NX-NH2 

161.1 

101.3 

pounds.  It  is  to  be  noted  that  experimental  determination  is 
often  unique  and  therefore  experimental  uncertainties  are 
often  underestimated.  The  same  method,  once  applied  to 
nitroso,  nitro,  nitrite,  and  nitrate  compounds  derived  from 
MMH  and  UDMIT,  leads  to  data  given  in  Table  2.  McQuaid 
et  al.  [12, 13]  used  another  ab  initio  method  to  compute  the 
gas-phase  standard  enthalpy  of  formation  of  some  of  the 
CH3-NX-NH?  species  given  in  Table  2.  Comparisons  of 
methods  are  always  interesting  especially  when  no  exper¬ 
imental  data  are  available.  Their  results  are  given  in  Table  3 
for  two  of  the  28  compounds  presented  in  Table  2.  Results 
given  in  Tables  2  and  3  arc  consistent  for  the  nitro  and  nitrite 
compounds  considered.  Therefore,  the  data  presented  for  26 
remaining  compounds  presented  in  Fable  2  can  be  consid¬ 
ered  to  be  reliable.  These  results,  and  the  reappraisal  of 
some  thermodynamic  data,  may  help  to  explain  hypergolic 
ignition.  All  the  compounds  in  Table  2  are  endothermic 
compounds,  i.e.  their  enthalpies  of  formation  arc  positive. 
Such  compounds  are  able  to  decompose  exothermically 
under  thermal  stress  and  this  may  be  of  interest  to  explain 
ignition  or  reactivity  characteristics  and  to  interpret  hard 
start  phenomenon  observed  occasionally  in  engines  using 
NTO/MMH. 

Another  theoretical  method  specifically  devoted  to  gen¬ 
eral  organic  compounds  not  containing  NO,  N02,  ONO, 
and/or  0N02  groups  was  then  used  to  compute  the 
thermochemistry  of  dimethyltetrazanes.  These  compounds 


Table  <1.  Gas  phase  standard  enthalpy  of  formation  at  298  K  of 
hydrazine,  MMH,  and  UDMH  calculated  by  using  a  ab  initio 
method  developed  for  this  study.  Units  are  kJ  mol-1. 


Compound 

Method 

Fxperimental 

used  here 

data 

Hydrazine 

95.4 

95.4 

MMH 

93.3 

94.1 

Table  5.  Gas  phase  standard  enthalpy  of  formation  at  298  K  of 
three  dimethyltetrazanes  by  using  a  ab  initio  method  developed 
for  this  study.  Units  are  kJ  mol-1. 


Compound 

Method 

Fxperimental 

used  here 

data 

CH3(NH2)NN(NH2)CH3 

269.4 

- 

CII3NHNHNHNHCH3 

CH3(NH,)NNHNHCH3 

258.6 

- 

262.3 

- 

are  formed  through  the  termolecular  recombination  of 
two  methylhydrazyl  radicals,  namely  CHiN*NH2  and 
CH3NHNH\  Their  termolecular  recombination  reactions 
lead  to  three  different  dimethyltetrazanes  according  to: 

CHjN'NHj  4  CH3N'NH2  4M- 

CH3(NH2)NN(NH2)CH3  4-  M 

CH3NHNH*  4  CH3NHNH*  4  M  - 

CH3NHNHNHNHCH3  +  M 

CH3NHNH*  4  CH3N*NH2  4  M  — ► 

CH3(NH2)NNHNHCH3  4  M 

The  agreement  between  experiments  and  calculations  is 
shown  in  Table  4  for  hydrazine  and  MMH.  Table  5  gives  the 
calculated  data  for  the  three  dimethyltetrazanes  given 
above.  The  formation  of  these  tetrazancs  is  exothermic 
and  these  compounds  are  endothermic,  i.e.  they  can 
decompose  exothermically.  T  hese  species  can  play  a  role 
in  ignition  and  this  role  needs  to  be  further  studied. 


4  Conclusion 

Protocols  have  been  derived  and  validated  for  predicting 
the  thermody  namic  data  needed  for  thermochemical  calcu- 
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la t ions  and  for  building  detailed  chemical  kinetic  models  to 
simulate  hypergolie  ignition  and  eombustion.  The  level  of 
theory  hybrid  density  functional  theory  with  a  modest  basis 
seh  plus  empirical  corrections,  is  found  appropriate  for  the 
current  state  of  knowledge  (quite  minimal)  of  their  ther¬ 
mochemistry,  This  paper  provides  Thermochemical  data  for 
31  molecules  potentially  formed  during  the  MMH/NTO 
reactivity.  No  experimental  data  are  available  for  these 
compounds.  The  model  previously  proposed  to  explain 
hypergolie  NTO/MMH  ignition  can  be  improved,  as  can 
almost  all  kinetic  models.  Detailed  kinetic  models  devel¬ 
oped  for  mixtures  easier  to  study  and  therefore  intensively 
studied,  such  asp  re  mixed  methane/airorhydrogen/airt  have 
been  continuously  improved  over  the  past  50  years,  and  new 
models  have  been  proposed  in  recent  years  to  explain 
ignition  characteristics  of  these  mixtures.  Therefore,  com¬ 
plementary  studies  need  to  be  performed  for  the  complex 
MMH/NTO  combustion  chemistry.  The  impact  of  the 
thermochemistry  proposed  in  this  study  on  the  predictions 
of  the  MMH/NTO  detailed  kinetic  model  must  he  exam¬ 
ined.  Alternative  reaction  paths  also  need  to  be  considered. 
The  development  of  visualization  techniques  and  of  new 
analytical  techniques,  along  with  the  development  and 
application  of  methods  for  calculating  the  thermochemistry 
of  reactants,  products,  and  intermediate  species  should  lead 
in  the  near  future  to  a  better  understanding  of  the  features  of 
hypergolie  ignition  and  to  new  tools  for  the  design  and 
screening  of  new  hypergolie  fuels. 
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